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ABSTRACT 

The motions of small moons through Saturn's rings provide excellent tests of radial 
migration models. In theory, torque exchange between these moons and ring particles 
leads to radial drift. We predict that moons with Hill radii r H ~ 2-24 km should 
migrate through the A ring in 1000 yr. In this size range, moons orbiting in an empty 
gap or in a full ring eventuall y migrate at the sam e rate. Smaller moons or moonlets 



such as the propellers (e.g., iTiscareno et al.ll2006l ) - are trapped by diffusion of disk 



material into corotating orbits, creating inertial drag. Larger moons - such as Pan 
or Atlas - do not migrate because of their own inertia. Fast migration of 2-24 km 
moons should eliminate intermediate-size bodies from the A ring and may be responsible 
for the observed large-radius cutoff of r H ~ 1-2 km in the size distribution of the 
A ring's propeller moonlets. Although the presence of Daphnis (r H ~ 5 km) inside 
the Keeler gap challenges this scenario, numerical simulations demonstrate that orbital 
resonances and stirring by distant, larger moons (e.g., Mimas) may be important factors. 
For Daphnis, stirring by distant moons seems the most promising mechanism to halt 
fast migration. Alternatively, Daphnis may be a recent addition to the ring that is 
settling into a low inclination orbit in ~ 10 3 yr prior to a phase of rapid migration. We 
provide predictions of observational constraints required to discriminate among possible 
scenarios for Daphnis. 

Subject headings: planetary systems - solar system: formation 



Introduction 



Migration is the steady radial drift of a massive planet or satellite through an astrophysical disk 
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planet lies embedded in a disk of gas or solid particles, it tries to clear material from its orbit. Clear- 
ing by a massive planet creates wakes inside and outside its orbit. These density perturbations exert 
a torque on the planet, which produces an inwards or outwards migration through the disk. Analytic 



timescale much shorter than the disk lifetime ( 


Goldreich & Tremaine 


1979; 


Lin &; PaDaloizou 


1979; 


Goldreich & Tremaine 


1980; 


Ida et al. 


2000; 


Tanaka et al. 
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Papaloizou et al. 


2007; 


Kirsh et al. 



200£ 



This process provides a popular explanation for the architectures of many planetary systems. 
Migr ation plausibly accounts for t he orbits of some "hot Jupiter s" very close to their parent stars 



(e.g., iMasset &; Papaloizou! 120031 ; lAlexander &; Armitage 1 120091 . and references therein) and the 



compact orbital configurations of several multiplanet sys tems, including Kepl er 11 which has five 
super-Earths orbiting within 0.3 AU of a Sun-like star (jLissauer et al.ll201ll ). In our own solar 
system, migration appears to be a key dynamica l mechanism for arranging the orbits of the outer 



planets (|Malhotra 



1993 



Hahn Malhotra 



1999 



Levison et al.1 120071 : iMorbidelli et al.ll2007h . 



Despite its broad applicability in exoplanetary systems, testing the predictions of migration 
theory is challenging. In a young planetary system, growing protopla nets gravitationally stir par - 
ticles in the disk and limit the formation of coherent density wakes (jBromlev &; Kenvonll2011bl ). 
For a single planet, m igration through a gaseous disk also depends on the dynamical and thermal 



state of the disk (e.g.. iPaardekooper et al.ll2010l . 120111). Real disks contain several planets in a con 



stant ly evolving mix of gas and solid particles (e.g., iBromlev &; Kenvonll2011al ; lYoudin Kenvon 
2012 , and references therein). Current calculations only partially address this complexity (e.g., 



Masset fc Snellgrove 



2001 



Morbidelli et al.ll2007l : lD'Angelo & Marzarill2012l ). Thus, theory is hard 



pressed to derive robust predictions applicable to observed systems. 

Saturn's rings provide an interesting laboratory in which to test migration theory. The A ring 
is a geometrically thin disk of solid particles with orbital periods sh orter than a terrestrial day. 
Aside from countless cm- to m-sized particles (e.g., ICuzzi et al.l 12009!) . the rings contai n many 0.1 
to 1 km moonlets that could easily migrate through the disk (ITiscareno et al.ll2010l ). NASA's 
Cassini-Huygens mi ssions have also revealed a variety of density perturbations within the rings 
(jColwell et al.ll2009r). identified a herd of "pro peller" moonlets (named for the wakes of ring par- 



ticles they creat e 



( Showalter 



1991 



Tiscareno et al 



Porco et al 



2006 . 120081 ) , and discovered the small moons Pan and Daphnis 



20051 ), Together, these features allow robust tests of our understand 



ing of gravity and dynamics in a complicated system of solid particles (see also ICuzzi et al. 
Espositdl2010l . and references therein). 



2010 



Previous theoretical analyses of migration in Saturn's rings have concentrated on the non- 
Keplerian, radial motion of propeller moonlets within the A ring (jTiscareno et al.ll2010l ). These 
moti ons may result from interactions between m o onlets and stoch astic density fluctuations in the 
disk (jCrida et al.1 bold : iRein fc PapaloizoJ bold : IPan et al.1 hold) or between moonlets and ring 



particles on loosely bound "frog" orbits around the moonlet (|Pan Chiang! |2010| . |2012| ), These 
phenomena, or perhaps direct collisions between moonlets and large ring particles, may cause the 
observed radial motion in some moonlets. On the other hand, large-scale r ing structure, pa rticularly 
radial density waves, may provide torques that suppress radial motion (jTiscarenol 120121 ) . In this 
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picture, the observed non-Keplerian orbits come from recent kicks received by moonlets, which then 
settle back to their original Keplerian configuration. Fortunately, each of these theoretical scenarios 
makes a unique prediction for the e volution of moon let orbits. Upcoming Cassini observations will 
be able to distinguish among them (jTiscarend 12012 ). 

Here, we focus on other modes of migration for moonlets and small moons embedded in Saturn's 
17,000 km wide A ringf^|. Aside from stochastic interactions, larger moonlets and small moons like 
Daphnis are candi dates for "fast" m igration modes that satellites experience in a dynamically cold, 
planetesimal disk (jlda et al.1 l2000l ). Because fast migration is a critical component of migration 
theory in exoplanetary systems, understanding fast modes in Saturn's rings improves our ability to 
predict migration in a broad range of environments. 

We begin with a theoretical overview of migration in a particle disk (§2) and then specialize 
the theoretical estimates to conditions appropriate for Saturn's A ring (§3). In §4, we confirm the 
analytic estimates with numerical simulations, but demonstrate that stochastic effects, resonances 
with more distant moons, and stirring by distant moons can slow down or halt migration of small 
satellites through the ring. In §5, we examine the migration of Daphnis, which is located within the 
Keeler gap. If Daphnis only interacts with ring particles, its migration should be fast and detectable 
with current satellites. However, orbital resonances (e.g., Prometheus 32:31) and stirring by moons 
out of the ring plane (e.g., Mimas and Tethys) may have a strong impact on Daphnis' orbit and 
migration rate. Measuring the evolution of Daphnis' orbit can constrain these possibilities. We 
conclude in §6 with a brief summary. 



2. Migration in a particle disk 



To test migration theory in Saturn's rings, we develop basic predictions for an idealized disk 
of solid particles orbiting a much more massive planet. For any satellite within the disk, the Hill 
radius, r H , defines a spherical surface where the gravity of the satellite roughly balances the gravity 
of the planet. This radius also characterizes the annular width of the satellite's corotation zone, 
an annulus where disk particles on nearly circular orbits oscillate back and forth in the satellite's 
rotating reference frame (e.g., on ho rseshoe orbits). Here we ta ke the halfwidth of the corotation 
zone to be X co r H , with X co = 2 (e.g.. iBromlev Kenvonll2011bl ). When the corotation zone is full 
of small particles, migration is not efficient. 



Two physical processes set the density of smaller particles in the corotation zone. The massive 
satellite gravitationally scatters smaller particles out of the zone. Viscous spreading drives the 
diffusion of smaller particles back into the zone. When scattering overcomes diffusion, the satellite 
opens up a physical gap in the disk. For convenience, we define r gap as the Hill radius of the 
smal lest satellite capable of clearing out its corotation zone (see also 



Lissauer et al 



1981 



Rafikov 



2001 



). When the corotation zone is nearly clear of smaller particles, migration can be very efficient. 



In the spirit of the definitions in iTiscarenol |2013h . we use "moonlet" to describe any of a numerous and perhaps 
unresolved population with physical radii of < 1 km and "small moon" to refer to the larger objects Daphnis, Pan or 
Altas. 
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When a satellite has r H > r gap , it undergoes fast migration. In this situation, the satellite can 
gravitationally pull corotating material all the way past its corotation zone, hauling itself radially 
inward or outward like a rope climbe r. In gaseous disks, this mode of migration is often called 
'type III' (jMasset fc Papaloizoul 120031 ) to distinguish it from the slow, type I migration of planets 
incapable of forming a gap a nd the intermediate pace of typ e II migration of more massive planets 
surrounded by a large gap@ (jWard 1997 ; Brvden et al. 1999 ). 



For more massive satellites with r H 3> ''gap; fast migration stalls. As the mass of a satellite 
grows, its Hill radius - and the width of its corotation zone - also grows. At some point, the satellite 
cannot scatter particles out of the corotation zone fast enough: some disk material remains inside 
the corotation zone and interacts with the planet. This interaction reduces the torque and slows 
migration. We define rf ast as the Hill radius where fast migration begins to stall. Thus, satellites 
with r gap < r H < rf ast undergo fast migration. 



In physical units, the three critical Hill radii for fast migration are: 

[Hill radius], 



gap 



''fast 



a UmJ ' 

0.4 (^aT) 1 / 3 , 
1.7a 



M 



[minimum ?* H to form a gap] , 
[maximum r H for fast migration] 



(1) 
(2) 
(3) 



In these expressions, m is the mass of an object orbiting with semimajor axis a and period T around 
a central star or planet with mass M. The disk has surface density X and radial viscosity f r ad- 

Deriving r gap and rf ast relies on identifying important timescales in the disk. Rates of migration 
are set by the satellite's Keplerian orbital period T. Although an oblate planet and gravitational 
perturbations by surrounding massive moons prevent a satellite from executing standard Keplerian 
orbits, the Keplerian orbital period is a useful fiducial for all other timescales. Close encounters 
between the satellite and particles occur every synodic period T syn , the timescale for particles on 
orbits with separation Aa to return to the same orbital phase. Massive satellites scatter particles 
away from their orbits, opening a gap of half-width Aa within the disk clearing time , T r \ par . For 
parti cle disks, the disk clearing time for a small gap is a few times larger than T syn (e.g. jKirsh et al 



20091 ) . Viscous processes between disk particles fill gaps on the viscous (or gap-filling) time 



which depends on the disk's radial viscosity i/ ra d (e.g.. lLissauer et al 



1981 



Hourigan &; Ward 



198 



Physical collisions damp the velocities of disk particles on the damping timescale Td amp . Col- 
lisional damping is a random walk- type process that depends on A ra d, the radial distance between 
collisions, and /i OS s> the fractional loss of kinetic energy per particle in a collision (equivalent to 
the s quare of the coefficient of restitution). The mean- free- path A is related to the optical depth r 



(e.g., iCook Franklin 



1964 



Goldreich fc Tremaindll978al ). In the low optical depth limit, a ring 



2 In a gaseous disk, uncertainties is the disk viscosity (and hence the radial and vertical gradients in the pressure, 
temperature, and other physical variables) make migration calculations very challenging. Thus, it is uncertain whether 
type III migr ation is pos s ible in a gaseous disk. In a particle di s k, tho u gh, Type III migra tion is generally accepted. 
Startin g with Malhotra ( 19931 ) and continuing with Ida et al. ( 200ol ) . Kirsh et al. ( 20091 ). and Bromlev fc Kenvonl 
( 2011 bh . fast/rapid/runaway migration is a standard phenomenon in particle disks. 
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of monodisperse particles with radius r p has r = 3S/4p p r p , where p p is the particle mass density 
(jGoldreich fc Tremaindll978al : ISalmon et al.ll2010l ). A more general, realistic treatment of optical 
depth requires consideration of the particle size distribut i on and coherent structures that may form 
as a result of self-gravity (e.g.. lSalo Karialainenll2003l ; iTiscarenol 120131 ) . 



For Aa <C a, these timescales are: 

.1/2 



a 3 

UM 



J- syn 
^clcar 

-^damp 



2?r 

2aT 
So 

AT 

^ l - L syn 
(2Aa) 2 

^rad 
A 2 . 



4r p p p a 



./"loss^rad ^JloSS^ 



Kcp 



[orbital period], 

[synodic period for orbital separation Aa], 
[gap clearing time (half- width Aa)], 
[gap filling/ viscous diffusion], 
[damping time]. 



(4) 
(5) 
(6) 
(7) 
(8) 



To derive the physical scales for fast migration, we compare the timescales for the relevant 
processes. When T c i ear = Tflii, a satellite clears its corotation zone as fast as viscous spreading 
fills the zone back up. This equality defines r gap , the smallest Hill radius for a satellite capable 
of opening a gap in the disk. Balancing the torque from the satellite on disk material with the 
torque from the disk on th e satellite is roughly equivalent to setting r r i P? r = Tm and yields similar 



Ward 


1997; 


' gait 
Rafikov 


\^->-, 
2001) 



Lissauer et al 



1981 



Hourigan fc Ward! 1 1984 : IWard fc Hourigan 



1989 



Aside from establishing which satellites can open gaps in the disk, r pa p s e ts a lower limit on the 



Hill radii of satellites capable of undergoing fast migration (jlda et al.ll2000l : iMasset h Papaloizou 



2003 



Kirsh et al.l 120091 : iBromley fc Kenyonl l2011bl ) . When a satellite pulls disk material near the 



edge of its corotation zone across its orbit, it recoils in the opposite direction. If the net recoil is 
large enough, the satellite does not interact with scattered particles when they pass by the satellite 
one synodic period later. Satellites with r H < r gap cannot clear their corotation zones and do not 
recoil enough to overcome the inertia of material along their orbits. Satellites with r H 3> r gap do 
not recoil enough to avoid a second encounter with scattered particles. 

To determine the range of Hill radii for satellites that undergo fast migration, we identify the 
maximum Hill radius where the time for a satellite to move a radial distance equal to the width of 
its corot ation zone is le s s than the synodic period of particles at the edge of the corotation zon e 



(see also llda et alJl200(l : Masset fc Papaloizoull2003l : iKirsh et al.ll2009l : iBromlev fc Kenvonll2011bT ). 
To derive this Hill radius^ which we define as rf ast , we set r syn = rir^ ^da^t/dt} -1 where da^t/dt 



is the fast migration ratqj (e.g., eq. 12 below; see also llda et al.ll2000l ). To avoid repeat encounters 
with material on horseshoe orbits within the corotation zone, we set r\ = 2. Although choosing r\ = 
3.5 would eliminate all possible repeat encounters, material trapped on horseshoe orbits provides 
most of the drag on a drifting moonlet. Thus, we choose 77 = 2. 



For disks of solid particles, the viscosity f ra d depends on kinematics. ^From lGoldreich fc Tremaine 



Bromlev fc KenvorJ (2011a) derive this rate from the recoil due to a single scattering event and an encounter 



frequency integrated over the width of the corotation zone. 
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(|1978al . and references therein), the viscosity is a function of the optical depth, the orbital period, 
and the radial velocity dispersion of particles v p : 



^rad 



4?r 1 + t 2 



(9) 



When r is small, i/ rac j oc r; when r is large, u Iai ^ oc r 



In addition to the surface density and the visc osity, the vertical scale height of the disk h also 
plays an important role in satelli t e migration (e.g.jLin &: Pap aloizou 1979; Gold reich fc Tremaine 



1980; Ward 



19971 : llda et al. 



2000 



Papaloizou et al.1 120071 : iBromley fe Kenyon 



201l"bfh For a. disk of 



particles, the scale height depends on the vertical velocity dispersion v z : 

h = v 7 Q~ l , 



(10) 



where £1 is the lo cal angular vel ocity of particles. If disk particles follow a Rayleigh distribution, 



Vp/2 (e.g., lOhtsuki I Il992l . and references therein). The local (3-D) yiscosity then scales 



approximately as y oc h 2 as in the a prescription for turbulent gaseous disks (jShakura fe Sunyaev 



19731 : iPringldl 198 ll ) . The true relationship between v z and the mean 3-D velocity dispersion is more 
complicated, espe cially in Saturn's A r ing where self-gravity wakes can increase velocities in the 
plane of the ring (|Tiscareno et al .1 120071 ). 



With these definitions and physical expressions for critical parameters, the theoretical radial 
drift rates for the three modes of migration are: 



dQcmb 
dt 

dafast 
dt 

dt 



M h 2 a T 

5.37ra 2 S a 
1 M T 

7ro 2 S r u a 
-16- 



r H < h, r gap [embedded, no gap; Type I] (11) 



h, 7-g a p < r H < r fas t [fast; Type III] 



(12) 



M a T 



3^rad 

2a 



r H ~ r gap 



[embedded+gap; Type II]. (13) 



When the Hill radius is smaller than the disk scale height and the gap radius, the satellite has a 
very weak gravity, cannot open up a gap, and produces very weak density perturbations in the disk. 
Only very sl ow, embedded migratio n is possible. Its rate (eq. [11]) follows from weak scat tering 
theor y (e.g.. lLin &; Papaloizou! 1 1973 ) and is analogous to Type I migration in gaseous disks (jWard 
19971 ). The scale height is important in setting this rate because it is a measure of the distance at 
which the ring pa rticles transition fro m random motion to coherent, shear flow from the satellite's 
perspective (e.g., iTanaka et al.1 l2002l ). F ast migration (ea . [T2| ) is the ideal case of a satellite 
spiraling through a dynamically cold disk (jMasset &; Papaloizou! 120031 ). The rate is independent of 
satellite mass; the greater gravitational reach of a more massive satellite compensates for its higher 
inertia. The migration rate for r H > r tast is the Type II analog for a planet which opens a gap 
in a particle disk (eq. p3]). In the low viscosity limit, differential torques arise from materi al at 
the edges of the gap, which lies at or beyond the corotation zone ([Bromley &: Kenvonll2011bl'). In 



the high viscosity limi t, the migration rate is set by the viscous timescale (see, e.g., IWard 
Papaloizou et al.l 120071 ). 



1997 
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The transiti on from type I migration with no gap to fast an d type II migration with a gap 
is complex (e.g., iHourigan &: Ward! 11984 : IWard &: Houriganl Il989l ). When a small satellite with 
r H < r gap accretes disk material and grows to larger radii, it drifts through the disk more rapidly (eq. 
[lip- If this satellite drifts across its corotation zone more rapidly than it can clear a gap, it drifts 
through the disk with a partially cleared corotation zone. Because the drift time scales more weakly 
with mass (oc m~ l ) than the time to clear a gap (oc m~ 2 ), there is a range of satellite masses where 
a sat ellite could clear a gap if it could stop drifting (jHourigan fc Ward! 11984 ; IWard k, Hourigan 
19891 ). This range of masses roughly corresponds to the range where fast migration operates. 



The transition from type I to type II migration also depends on the scale height and eccentricity 
e of disk particles. When disk particles have large e, the drift rate slows by a factor ~ (r H /e p a) 3 
relative to the fast rate (see eq. [12] : Bromley & Kenyon 2011b). When h < r H < r gap , simulations 
suggest that nonlinearities and other competing effects - such as the inertial drag of material stuck 
in the corotation zone - substantially reduce the migration in this "transgap" range of Hill radii 



(ID'Angelo et al 



ring (ITiscareno 



2003). This conclusion is strengthened by a recent analysis of moonlets in the A 
20121 ) . which shows how radial density structures in the ring can produce torques 



that damp radial drift. 



For larger satellites with r H ~ ?"f ast , simulations ([Bromley &; Kenvonll2011bh also suggest that 
fast migration falls off at least as rapidly as (r^t/r^) 3 (see also lSver Clarkdll995l ). As r H increases 
(r H > r gap ,rf ast ), the satellite can maintain a gap in the disk and can exp erience migration with a 



fully- formed gap (eq. [13] ) . In the inviscid limit of this type II migration (jWardl 119971 ; ICrida et al 



20061 ). the rate is given by weak scattering theory. In a sufficiently viscous di sk, material piles u p 
at the gap edges. Then, the timescale for radial drift is the viscous timescale (IBryden et al.lll999l ). 



The sketch in Figure [T] summarizes the various migration modes. In addition to showing the 
relative magnitude of the rates in equations (|1 1H13|) as a function of r H , the sketch illustrates 
the impact of viscosity in the transgap regime and the high- and low-mass fall off of the fast 
migration mode. Other motions affecting a satellite's radial drift include stochastic migration (see 



Rein &: Papaloizou 



2010 



Crida et al. 



2010 



Pan et al.1 12012 ) and the pe riodic orbit perturba tions 



caused by material on corotating "frog" orbits centered on the satellite (jPan &; Chiangil2010l ). On 
yearlong timescales, these drifts (random walks with effective da/dt values that lie in the shaded 
region in the figure) can be comparable to linear-in-time radial migration in Saturn's rings. 



3. Migration in the A ring 



Saturn's A ring is a dynamically cold, thin disk of icy, meter-sized particles, lying in an annulus 
around Saturn that extends from a radius of 117,000 km to the Roche Division at 137,000 km. In 
our calculations, we adopt an o rbital radius of a = 130, 000 km as represen tative. For other physical 



parameters, observations (e.g.. ITiscareno et al. 



2007 



Colwell et al.l 120091 ) suggest: 



£ = 40 g/cm 2 
h = 10 m 

^rad = 85 cm 2 /s. 



(14) 
(15) 
(16) 
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Adopting these values allows us to specialize the results of ^21 to Saturn's A ring. 

In the A ring, the radii that delineate different modes of migration in equations ([TH3|) become 



1.7 



p 


1/3 - 


1 g/cm 3 





1.3 xlO 5 km 



^phys > 



gap 



Hast 



^rad 



24 



85 cm 2 /s 



a 



1/3 



a 



1.3 xlO 5 km 



1.3 xlO 5 km 



40 g/cm 2 



5/6 



km, 



1/2 



km. 



(17) 
(18) 
(19) 



With mean densities of ~ 0.5-1 g/cm 3 , the Hill radius and the physical radius of a moonlet or 
small moon in the A ring are similar. Prom equation (|18|) . satellites with Hill radii below a few km 
cannot open a full gap in the A ring. These bodies can only migrate in the slow, embedded mode. 
Satellites with r H > 20 km are too massive for fast migration. Between a few km and 20 km, small 
moons can migrate in fast mode. 



The timescales in equations (J3HSJ) specialize to 

1 3/2 



^clear 
Ti damp 



0.0015 



130 



a 



1.3 xlO 5 km 



a 



520 
3.7 
0.008 



1 km 




Aa 




1 km 




Aa 




" Aa " 


2 - 


1 km 





1.3 xlO 5 km 
a 

1.3 xlO 5 km 
85 cm 2 /s 



5/2 
5/2 



yr, 



yr, 



0.1 



loss 



^rad 
P 

1 g/cm 3 



yr, 

40 g/cm 2 



1 m 



3/2 



yr 



^ 1.3 xlO 5 km_ 

In the expression for the damping time, we adopt f\ oss = 0.1 (see ^4.6|) . 

Finally, the migration rates in the A ring, corresponding to equations (jllH13p . are: 



(20) 
(21) 
(22) 
(23) 
(24) 



dt 

^Qfast 

dt 
dt 



-10" 



±17 



4 


"10 m" 


2 


a 


1/2 


s 




h 




1.3 x10 s km 




40 g/cm 2 



km/yr 



a 



1.3 xlO 5 km 



-0.009 



24 km 



-i3/2 


S 




40 g/cm 2 




a 


1/2 


1.3 xlO 5 km 





km/yr 
S 

40 g/cm 2 



km/yr 



r H < 2 km 



2 km<r H <24 km 



r„ £ 24 km . 



(25) 
(26) 
(27) 



Equation (|25p gives the maximum embedded rate for r H ~ h, with the expectation that embedded 



satel lites with larger radii cannot migrate faster as a result of viscous effects (e.g., iD'Angelo et al 



20031 ) . Equation (|27|) expresses the migration rate in the limit of small viscosity. In a viscous disk, 



da v 



dt 



-3 x 10" 6 



1.3 xlO 5 km 



^rad 

85 cm 2 /s 



km/yr r H > 24 km, f rad > . (28) 
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These equations provide an idealized description of migration in the A ring. At small r H , we take 
the embedded migration rate for r H ~ 3h as a maximum value, since larger Hill radii are in the 
transgap regime where migration is likely suppressed. The large- r H , viscosity-dominated, Type II 
migration rate is formally below 1 cm/yr. 

To put these results in perspective, we estimate the timescales for migration through a distance 
of 10,000 km, approximately half the annular width of the A ring: 



Temb ~ 100 



1 km 


3 


" h ' 


2 


"1.3 x10 s km" 


1/2 


40 g/cm 2 


r H _ 




10 m_ 




a 







Tfast 



2,ap 



0.0006 



1.3 xlO 5 km 



-|3/2 



1.1 



370 



24 km 



a 

1.3 xlO 5 km 



40 g/cm 2 



Myr 



a 



1/2 



40 g/cm 2 



Myr 



85 cm 2 /s 

^rad 



Myr 



Myr r H < 2 km (29) 



2 km < r H < 24 km (30) 



r H > 24 km . (31) 



(viscous lifetime) . (32) 



Figure [2] summarizes results for the migration rate in the A ring as a function of Hill radius. 
Satellites in or near the A ring are included in the diagram. 

These timescales lead to several clear conclusions regarding migration through the A ring. 
Objects with Hill radii less t han roughly 0.25 km migrat e on a timescale comparable to the age 
of the solar system (see also iGoldreich &: Tremaind 11982 . and references therein) . This radius is 
roughly the maximum Hill radius of propeller moonlets in the A ring. Because the migration time 
scales as r^, fast migration may be responsible for the lack of propeller moonlets with r H > 1 km. 
Two moons, Pan and Atlas, have Hill radii very close to the fast migration cut-off rf ast . Although 
the predicted lifetimes for these moons are short compared to the age of the solar system, they are 
very long compared to the fast migration lifetime. 

The small moon Daplmis clearly provides a more extreme challenge to migration theory. Its 
radius lies right in the middle of the expected range of radii for fast migration. With an expected 
migration time of less than 1000 yr, either (i) Daphnis is a recent addition to the rings, (ii) the 
migration theory is incorrect, or (hi) additional factors contribute to its apparent lack of migration. 
In the following sections, we address these possibilities and their impact on the rates shown in 
Figure El 



4. Migration of small moons: N-body experiments 



To explore the migration of small moons through Saturn's A ring in more detail, we examine 
a set of A-body calculations. To perform these calculations, we use Orchestra, a hybrid A-body- 
coagulation code which includes additional module s for the radial diffus i on of gas and solid particles 
throu g h a d i sk sur round ing a star or a plane t . In IBromley &: Kenyonl (|2006l ). iKenyon &: Bromley 
(|2004l . 120081 . [2009j), and IBromley &: Kenyonl (|2011al ). we describe each component of Orchestra, 
including comprehensive results with standard test problems. Here, we provide a short overview, 
with a focus on the A-body component, which we use exclusively in the calculations for this paper. 
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4.1. N-body simulation method 

The iV-body algorithm distinguishes among three types of particles: the massive iV-bodies, 
the "swarm" (superparticles) , and the tracers. The massive iV-bodies gravitationally interact with 
all particles in a simulation. The swarm particles interact with and can influence the motions of 
the massive iV-bodies but do not interact with each other or with the tracers. The tracers respond 
to the gravity of the massive iV-bodies. Here, we use massive iV-bodies to track the motions of 
Saturn, its major moons, and the smaller moons. Unless noted otherwise, we represent the gravity 
of the ring using swarm particles. Thus, the ring is not self-gravitating. 

Orchestra uses a 6 th -order accurate symplectic integrator for all massive bodies in a compu- 
tational domain that encompasses the full 3-D extent of all orbits. In the case of the relatively 
numerous tracer or swarm particles, a software switch allows for 4 th -order accurate integration over 
a timestep that is typically l/300 th of the orbital period. During this step, the orbits of massive 
particles are interpolated in a way that allows them to respond to the gravity from swarm parti- 
cles. The code handles close encounters with a massive body by solving Kepler's equations. These 
features allow us to track up to 10 s ring particles for thousands of orbits, or 10 5 ring particles for 
millions of orbits. 

Orchestra can also accommodate rates of change in orbital elements, given by da/dt (semimajor 
axis), de/dt (eccentricity), and di/dt (inclination), from physical interactions among the swarm 
particles. Here, we use this feature to include the evolution of solid ring particles in a viscous disk. 
To calculate changes in these elements over a time step At, we use 



Aa = sgn(n — 0.5) VtVadAt (33) 
Ae = -e-^- (34) 

damp 

Ai = (35) 

damp 

where u is a standard uniform (pseudo)random variate used to implement a random walk to mimic 
the radial diffusion of the swarm particles. For individual isolated particles scattered within roughly 
a Hill radius of a small moon, the code has a switch to turn off this feature. 

In the calculations for this paper, we chose a damping time of 0.01 yr, comparable to the 
predicted timescale in equation (I24|) . The diffusion of small ring particles is an important aspect 
of the interactions between massive iV-bodies and the swarm. As long as it is within an order of 
magnitude of our choice, the precise value of Td am p does not make any noticeable difference in the 
outcomes of our simulations, 

To test this aspect of our iV-body code, we considered the evolution of a narrow ring of 
small parti cles evolving solely b y diffusion and collisional damping. Starting with initial conditions 



identical to I Salmon et al.l (|2010l ) , our simulations yield a gradual expansion of the narrow ring into 



a broad disk. For at least 100 Myr, the time e volution of the surface density distribution in our 



simualtions closely follows ISalmon et al.l (|2010l ). who solve the ID radial diffusion equation for a 
disk of small particles. Comparisons with our own solutions to the radial diffusion equation further 
demonstrate that the iV-body code treats particle diffusion and collisional damping accurately. 
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Finally, we developed and tested algorithms to model Saturn's oblateness, which can signifi- 
cantly affect orbital phases and inferred orbital elements for all satellites and ring particles. When 
we include Saturn's oblateness in the calculations, we set the gravitational acceleration at a distance 
x near the ring to 



d 2 x GM 



dt 2 x 3 



R 

l + 3J 2 P2(sin(/>)^ 



(36) 



where 4> is the latitude relative to the ring, P2 the second-order Legendre polynomi al, and the grav- 



itatio nal constant J2 = 0.01629 for an adopted radius of Saturn R$ = 60, 330 km (jjacobson et al 



20061 ). 



For this paper, we consider two sets of simulations. When examining a simple physical process 
- such as a small moon migrating through a pristine ring - we treat the dynamical interactions 
between a small moon and ring particles using swarm particles and the damping/viscous interac- 
tions among ring particles (eq. [35]), but we do not generally include the oblateness of Saturn or 
interactions with distant moons (e.g., Mimas, Rhea, Dione, etc). Eliminating oblateness and the 
distant moons allows us to investigate a larger range of input parameters per unit cpu time. Once 
we establish a typical behavior, we verify that including oblateness and interactions with distant 
moons do not change this behavior. In more complex simulations (e.g., for Daphnis in §5), we 
include oblateness and the interactions with distant moons and note the included moons in the 
text. 



4.2. Comparison with other techniques 



The iV-body algorithm in Orchestra differs from 'patch' simulations of Saturn's rings, where 
the computation al domain is a comoving rectilinear patch centered on a circular Keplerian orbit 
in the ring (e . 



Perrine et al. 



Lewis &: Stewart 



2011 



2000 



Porco et al. 



2008 



Salo &: Schmidt 



2010 



Lewis et al 



2011 



Patch calculations include 10 5 to 10 6 ring particles with typical masses of 
5 — 10 x 10~ 5 relative to the satellite of interest. For patches with sizes of roughly a square 
kilometer, these calculations can track the interactions and long-term evolution of ring particles 
with realistic sizes, ~ 1-10 m. To understand the origin of local structures within the rings, these 
calculations are ideal. However, periodic boundary conditions at the 'ends' of a patch limit the 
ability of these calculations to include perturbations from material outside the patch. Thus, it is 
challenging for these simulations to track the response of a moonlet or small moon to the global 
evolution of ring particles. 

In Orchestra, the iV-body code follows much larger particles, ~ 10-100 m, over complete orbits 
around the central mass. As long as the satellite mass exceeds the ring particle mass by a factor 
of a few hundred, calculations suggest that these superparticles can accuratel y represent smalle r 
ring particles and their effect on an embedded satellite (Fig. [3] below; see also iKirsh et al.1 120091 ) . 
This approach also enables straightforward inclusion of resonances and other perturbations from 
distant moons. In our simulations, the number of superparticles is too small to identify small-scale 
or short-term structures within the ring particles. Thus, these calculations accurately follow the 
evolution of a small moon at the expense of a large graininess in the structure of the ring. Because 
our focus is on the long-term evolution of a satellite in response to the entire ring and distant 
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moons, this compromise is reasonable. 



Although most previous TV-body work has focused exclusively on ring particle dynamics, sev 



eral studies conside r radial drift of satellites within or outside the ring (jLewis &: Stewart 



2009 



Charnoz et alJl201ll ; lTiscarendl2012l ). Despite their success in id entifying sources of ra dial drift, 
none of these simulations can track fast migration. For example, ICharnoz et al.l (120 111 ) elegantly 
track Type II-like migration by summing the torque on a moonlet at Lindblad resonances. However, 
the relevant p art of the ring for fast migr ation is closer to the corotation zone than the Lindblad 
resonance (see iMasset k, Papaloizoull2003l . for the analogous case of Type III migration in a gaseous 
disk). Thus, this approach will never yield fast migration. 



Lewis fc Stewartl (|2009l ) provide an excellent demonstration of chaotic drift of a moonlet by 



collisions and interactions with massive particles and self-gravitating wakes in the rings. However, 
their patch-like computational domain spans only a small fraction of a moonlet 's orbit around the 
planet. To produce fast migration in this type of calculation, we speculate that it is necessary to 
account for the phase and timing of a ri ng particle's orbit when it reenters the patch. In all of their 
calculations, the moonlets of interest to iLewis &: Stewart! (120091 ) are too small to open the annular 
gap in the ring necessary for fast migration (e.g., r H < r gap ). Thus, fast migration is not expected 
in any of their simulations. 



4.3. Onset of fast migration 

Sustained fast migration of a small moon requires one simple condition. The moon must be 
able to draw/scatter ring material from one edge of its corotation zone to the other and then drift 
a sufficiently large radial distance in the opposite direction to avoid encountering the scattered 
particles again. Moons that are too large cannot drift far enough before re-encountering scattered 
particles. Very small moons cannot clear enough of an annular gap and are inertially pinned by 
other corotating material at smaller orbital separation. Thus, moons which are too small or too 
large migrate very slowly. 

Despite this clear picture for sustaining fast migration, initiating fast radial drift is more 
subtle. The onset of fast migration requires (i) a satellite with a fairly empty corotation zone 
and (ii) enough disk material close by, within ~ 3 — 4r H , where dynamical interactions can draw 
material into the satellite's corotation zonjj|. Among plausible starting points, we consider two 
extreme initial conditions, a satellite lying in a completely empty gap in the disk and a satellite 
embedded in a continuous disk. 

Although a satellite within a continuous disk has enough material nearby to begin fast mi- 
gration, it must first clear its corotation zone. For satellites with r H > r gap> the timescale for 
scattering particles out of the corotation zone is smaller than the timescale for viscous diffusion to 
bring particles back into the corotation zone. Thus, the satellite begins to clear its corotation zone 



4 In addition to particles on horseshoe orbits, m aterial on initia lly circular orbits with Aa as large as 2\/3V H can 
be deflected onto crossing orbits with the satellite ( Gladman 19931 ). 
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(see also iLissauer et al.lll98ll ; iHourigan fc Wardll984 ; lRafikovll200ll ). Whether or not the satellite 
drifts slowly in the type I mode, any imbalance in the net amount of scattering into the inner 
disk relative to the outer disk creates a net drift sufficient to begin fast migration. As long as the 



physical properties of the disk and the satellite maintain the conditions T c i ear 
the satellite continues fast migration. 



> T m and r H < 



'"fast ) 



For a satellite within a large gap, the onset of fast migration is more delicate. Despite its 
empty corotation zone, the satellite must reach across the gap and draw material from the gap 
edges into the region, Aa < 2\/3r H ~ 3.5r H , where orbit crossings can begin. For a satellite with 
r H S> rfast) the half-width of the gap is much larger than ~ 3.5 r H . Disk material is too far away; 
the satellite can only migrate in the type II mode. For satellites with r gap < r H < rf ast , the edges 
of the gap are within ~ 3.5 r H . Material can then be drawn into crossing orbits with the satellite. 

Although satellites with r gap < r H < rf ast have fairly small gaps, it is possible to construct a 
disk where orbit crossings never occur. We first consider a disk where the surface density at the 
inner edge of the gap is comparable to the surface density at the outer edge. When (i) a disk has 
a vertical scale height h = 0, (ii) the satellite lies on a perfectly circular orbit centered in a gap 
with half-width Aa > 4r H , and (iii) there ar e no external torques, th e torque balance between the 
satellite and the gap edges is perfect (e.g., IHourigan &: Ward! 11984 ) . The satellite cannot draw 
material into its corotation zone and can only migrate slowly in the type II mode. 

In a real disk, the finit e scale height breaks th e symmetry, allowing the satellite to draw 
material from the gap edges (IHourigan &: Wardlll984l ). Once the satellite begins to draw material 
from the edges of the gap across its corotation zone, additional jostling of the gap edges or the 
satellite's orbit - for example, by density fluctuations in the disk or by gravity fluctuations due to 
distant moons - places the satellite closer to material at one edge of the gap than the other edge. 
The satellite draws more material into its corotation zone. Fast migration commences. 

To try to maintain the satellite on a stable orbit, we can specify the surface density on the 
inner /outer edges of the gap to balance the slightly different angular momenta of particles on either 
edge of the gap. We can also specify that the number of orbit crossings from material in the inner 
disk precisely balance the number from the outer disk. Every ensemble of scattering events then 
leads to no radial motion. Neglecting internal ring dynamics and other external processes, the 
satellite can maintain a wide, clear gap and remain on its initial orbit indefinitely. 

In a real ring, however, this balance is metastable. Radial diffusion tries to fill the gap. 
Perturbations from scattering by self-gravitating wakes, direct collisions, or asymmetries in the 
Type I- like torques ( eq. |llj ) from more distant ring material pushes ring particles into the gap 
(e.g., iTiscarend 120121 ) . Jostling by distant moons also tends to fill the gap. If these processes fill 
both sides of the gap identically, it is possible to maintain the balance. If these perturbations leave 
a few extra particles on one edge of the gap, the additional crossing events will pull the moon in 
that direction. Once the moon is drawn to either gap edge, it will pull even more particles from 
that edge, enhancing its radial drift. Fast migration again commences. 



To explore these two possibilities in more detail, we examine a suite of numerical simulations 
for satellites in disks with physical characteristics similar to Saturn's A ring. We begin with a 
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discussion of fast migration for a satellite in a completely empty gap. We then consider a satellite 
in a continuous disk. 



4-3.1. Fast migration in smooth rings 

To illustrate metastability of a satellite in a pristine gap, we set up an initially stable configu- 
ration with balanced torques from either edge of the gap. The simulations consider two plausible 
initial conditions for particles spaced uniformly in radius through the ring: particles spaced (i) ran- 
domly or (ii) uniformly in the azimuthal direction. The calculations also consider two types of ring 
particles: (i) massive swarm particles, where we calculate the radial drift from individual iV-body 
interactions directly or (ii) massless tracer particles, where we derive the drift from the change 
in the total angular momentum of the tracers. This approach allows us to test the sensitivity of 
the migration rate to the initial orbits of ring particles and to the method for calculating torques 
between the moon and ring particles. 

For a small moon with r H = 7 km and a gap with half-width = 2 r H , fast migration initiates 
in every simulation in this suite of tests (Fig. [3]). Moons interacting with tracer particles tend 
to begin migrating earlier than moons surrounded by massive swarm particles. Moons within a 
random distribution of ring particles appear to start migrating earlier than moons in a uniform 
distribution. In all cases, the onset of fast migration requires a short time, smaller than the synodic 
period of material at the gap edge. When the mass ratio of the moon to a ring particle is large, 
> 10 3 , the onset is slow and regular. When the mass ratio is < 10 2 , migration is more stochastic 
and start more quickly due to individual scattering events that push the moon toward the edge of 
the gap. Still, both sets of initial conditions, both approaches to calculate the radial drift, and all 
mass ratios yield roughly the same migration rate. 

To explore the onset of fast migration in more detail, we conducted additional tests with 
larger gaps and different diffusion timescales. The range of gap sizes, with half-widths Aa from 
2r H to 6r H , includes two nominal half- widths that are important from gap- formation theory (see 
the Appendix) : the first is Aa « 3r H , derived by balancing torques between the satellite and the 



disk (jLissauer et al.lll98ll ). while the second is Aa ~ 2.5r H , derived by balancing scattering with 
viscous diffusion. In all these tests, ring particles slowly diffuse into the gap. The width of the 
gap around the moon gradually decreases. With the nominal diffusion time of 10 yr for a small, 
gap-clearing moon, larger gaps fill in more slowy, and imply longer times before migration can 
start. When the gap width closes to roughly twice the corotation half-width, Aa ~ 4r H , the moon 
starts to chaotically scatter ring particles. Fast migration then ensues. With no diffusion, large 
gaps remain pristine. Although slow, type I torques can cause radial drift on very long timescales, 
fast migration cannot begin until the moon is within a corotation radius of either edge of the gap. 
Therefore, eliminating particle diffusion in a large gap effectively eliminates fast migration. 



- 15 - 



4-3.2. Fast migration in wide gaps: encroaching ring edges 

To provide more details on how a moon within a large gap begins fast migration, we consider 
another example. The simulation begins with a small moon, r H = 4.9 km, in a perfectly circular 
orbit within a large gap of half-width Aa = 20 km ~ 4r H . The gap is centered within a larger ring 
of particles with properties appropriate for the A ring, equations (j!4H16j) . As in other calculations, 
ring particles have a diffusion time of 10 yr. The simulation spanned more than 50 years, or 
25,000 orbits. A sequence of snap shots shows how gap edges are sculpted and how particle orbits 
diffuse close enough to the satellite to trigger fast migration. 

Figure S] illustrates how ring material slowly diffuses into a wide gap and eventually leaks onto 
orbits that interact strongly with the small moon. Before the moon starts to migrate, very little of 
this material enters the moon's corotation zone. Instead, slight imbalances draw the moon slightly 
closer to one gap edge than the other. This process is slow. In the top two panels of Figure [H the 
time difference is roughly 40 yr. Only slight changes are visible. The gap is still well defined. The 
wakes created by the moon are clear. 

Once changes begin, they are rapid. In Figure HJ the lower two panels show the scene only 
a few years later. Orbit crossing is rampant. As the moon drives into one side of the gap, it 
pulls material across, opening a new gap ahead of it and fillin g in the space behind it. T his new, 



smaller and less well-defined gap comoves with the satellite (see lBromlev Kenvonll2011bl . and the 
associated animation in a much broader disk). Although the original gap remains behind, viscous 
diffusion slowly fills it in. 

Tests with moons of other sizes in gaps of varying widths yield similar results. Diffusion tries 
to fill the gap; scattering tries to empty it. Eventually, there is an imbalance which pushes the 
moon towards one edge of the gap. The accelerating imbalance between scattering and viscous 
diffusion allows the moon to draw more and more material from the gap and to move closer and 
closer to the edge of the gap. Once the imbalance reaches a critical level, fast migration begins. 

The starting conditions for the results shown in Figure |4] are almost a realistic representation 
of Daphnis, the A ring and the Keeler gap. The single difference is the initial orbit of the moon. In 
the simulation, the initial orbit has e = i = 0. The real orbit is slightly eccentric and inclined to the 
ring plane. As we show in §5, the time required for a moon to begin fast migration depends on the 
initial e and i. However, the overall outcome is identical: the moon always begins fast migration. 



4-3.3. In situ formation and clearing the corotation zone 



When it grows by coagulation/accretion, a satellite forms within a dense disk of particles. 
With a full corot ation zone, this sa tellite cannot migrate in the fast mode. Previous numerical 



experiments (e.g., iKirsh et al.ll2009i ) demonstrate that a satellite can gravitationally scatter away 



enough of the corotating material to trigger fast migration in a particle disk. Here, we illustrate 
the same phenomenon for a small moon in Saturn's A ring. Figure [5] shows results for a simulated 
moon with r H = 14 km in the A ring. The figure tracks the orbital distance and the amount of 
material in the corotation zone as a function of time. After the mass in the corotation zone falls 
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to about a third of its starting value, fast migration starts. The mass of corotating material then 
peaks up as the moon plows through fresh material during its inward spiral through the disk. 

Figure [5] shows th at a small moon rapidly clears it s corotation zone (for other examples of this 
phenomenon, see also lRafikovll200ll : iKirsh et al.ll2009l ). Within 10-20 yr, gravitational scattering 
removes roughly half the mass of the corotation zone. During this period, the moon slowly moves 
radially inward and outward in response to perturbations from swarm particles in the corotation 
zone. After ~ 100 yr, the moon has removed enough material from the corotation zone to begin fast 
migration. After migrating inward roughly 200 km, the moon's progress is halted as it encounters 
the edge of the disk. 

Our simulations demonstrate that this behavior is independent of the general background 
perturbations from other small moons and larger distant ones orbiting Saturn. Comparison of 
calculations without distant moons (Fig. [5j black curve) and with a full complemento of moons 
(Fig. [5j grey curve) shows modest differences resulting from stochastic variations in the orbits of 
swarm particles but no large variation in the timing or rate of migration. Results for a large suite 
of simulations yield the same conclusion: fast migration of satellites with r gap < r H < rf ast is a 
robust phenomenon in the A ri ng. Thus, the general lack o f small moons in the ring is consistent 
with migration theory (see also iGoldreich &; Tremaindll982l ). 



4-3.4- Orbital settling into a ring 

The onset of fast migration depends on the initial orbital parameters of a satellite. When a 
small moon is scattered into an otherwise unperturbed region of a ring, it will probably have a large 
orbital eccentricity or inclination relative to the ring particles. The moon then tries to clear a region 
much larger than its corotation zone. When e > 10 -5 , the radial excursion of the moon exceeds 
the radius of its corotation zone. When % > 10" 7 , the moon spends most of its time outside the 
ring and cannot completely clear out mat erial along its orbit. However, gravitational interactions 



with the ring par t icles damp e and i (e.g.. IGoldreich &; Tremaine 



1988 



Hahn 



1981 



Borderies et al. 1984; Ward 



2007 



20081 ). Once e and i are small, the moon can clear its corotation zone and begin 



fast migration. 

To derive the damping time for a small moon with an inclined orbit, we perform a set of 
calculations with various initial e and i relative to swarm particles in the A ring. Figure [6] illustrates 
the rapid damping of e and i experienced by a small moon in a typical calculation. For moons 
capable of undergoing fast migration, the damping timescale is well within 1-2 decades. The 
inclination damps smoothly and reaches i < 10~ 6 in 10-20 yr. After moons reach e < 10 , the 
drop in e is more gradual and includes additional stochastic variations due to perturbations in the 
surface density of swarm particles. 



Once a small moon achieves a nearly circular orbit, the timescale to damp the orbit further is 



5 In these simulations, a 'full complement' of moons includes the ring moons Pan, Daphnis, Atlas, Prometheus, 
Pandora, Epimetheus, and Janus and the larger, more distant moons Mimas, Enceladus, Tethys, Dione, Rhea, Titan, 
and Iapetus. 
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much longer. When the inclination satisfies r H < ai, the moon can begin to clear its corotation zone. 
The timescale to reach this inclination is roughly the viscous time to fill the gap formed when the 
moon had much larger e and i. Figure [7] provides an illustration. A small moon, with r H = 6.5 km, 
experiences a slow damping followed by a transition to fast migration. 

Simulations of the orbital evolution of small moons with r gap < r H < rf ast and a broad range 
of starting e and i yield similar results. As the moons try to clear a gap much larger than their 
corotation zones, their e and i damp. While e and i damp, viscous spreading fills in the gap outside 
a moon's corotation zone. After i reaches a threshold, fast migration begins. 

The vertical scale height of the ring sets the threshold for fast migration of a small moon 
with an inclined orbit relative to the ring. The condition for fast migration is if as t < r H /a. When 
i > r H /a, moons cannot completely clear their corotation zones before viscous spreading fills them 
in. After damping reduces i to the threshold if as t, fast migration begins. 

We conclude that fast migration may be a common, if not ubiquitous, process in a disk similar 
in structure to Saturn's A ring. The mig ration time is extremely short compared to the ~ 0.1-1 Gyr 



lifetime of the ring (e.g.. lEspositoll2010i ). Thus, (i) all moons with r H > r gap are recent additions to 
the ring or (ii) physical structures within the ring prevent migration. To address these possibilities, 
we next consider how density perturbations within the ring impact migration. In §5, we examine 
the possibility that the moons are young. 



4.4. Fast versus stochastic migration 



The A ring is dynamically active, with a variety of long-li ved and transient p hysical structures. 
Man y features are as sociated with spiral density waves (e.g., ICuzzi et al.lll98ll ) or bending waves 
(e.g., IShu et al.lll983l ). In addition to the surface density enhancements immediately surrounding 



moonlets (the propelle rs) , nearby moonlets or small moons generate gravitational wakes in the rings 



wakes ( Hedman et al 



(IShowalter et al.lll986 ). Ad ditional variations of optical depth in the A ring result from self-gravity 



20071 ). which vary in structure across orbital resonances with Saturn's moons. 



Although simulating the full range of structure present in the A ring is beyond the current 
capabilities of Orchestra, we can establish how a variety of density perturbations impact fast mi- 
gration in the A ring. We begin with simulations designed to illustrate how dumpiness in the ring 
impacts migration. 

Figure [8] shows a suite of simulations of a moon with r H = 14 km at an orbital distance of 
130,000 km embedded in rings with a fixed surface density and with swarm particles of various 
sizes. The swarm particles in each simulation all have the same mass, but from run to run, that 
mass ranges from 0.2% to 10% of the moon mass. In all simulations, the corotation zone is initially 
free of disk particles in an annulus of 3r H . Fast migration begins nearly immediately. Without this 
initial clear ing of the corotati on zone, there is a time lag ~ T c i ear in the onset of fast migration (see 
Fig. 5 and kirsh et alibood ). 



These simulations demonstrate that density perturbations can change migration rates substan- 
tially. As the mass of a swarm particle increases, migration through the ring becomes more and 
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more erratic (Fig. [5]). When swarm particles are 10% of the satellite mass, migrati on begins to 



resem ble a random walk process, as envisioned in the stochastic migration scenario of ICrida et al 



(|2010l . see also Tiscareno 2012). 



For a more detailed illustration of discreteness effects in migration, we consider the migration 
of a small moonlet with r H = 700 m. As in Figure O the calculations include all of Saturn's major 
moons (e.g., Titan and Rhea) and smaller moons near the A ring (e.g., Prometheus and Pandora). 
To enable computational feasibility, we reduce the number of swarm particles by decreasing the 
surface density of the A ring by a factor of 100, setting the mass of swarm particles at 1-2% of the 
moonlet mass, and clearing the corotation zone of swarm particles. We also turn off viscous diffu- 
sion. Without viscosity, particles with Hill radii comparable to the propeller moonlets (r H < 1 km) 
can maintain a fairly empty corotation zone once they start migrating. 

Figure [9] shows typical results of these simulations. The noisy horizontal curve in light grey 
measures the fluctuations in semimajor axis - as measured in a reference frame centered on Saturn 
- of an isolated moonlet orbiting the planet but not migrating through the disk. These fluctuations 
result from the constant gravitational jostling from distant and nearby moons. When the moonlet 
is allowed to migrate (as in the black and dark grey curves, which have particles that are 1% 
and 2% of the moonlet mass, respectively), its semimajor axis shifts in a stepwise fashion which 
tracks the theoretical prediction (dashed curve). Over a one year time frame, the total migration 
is comparable to the variation in a from jostling by distant moons. 

These two simulations illustrate the difficulty of measuring migration rates in Saturn's rings 
and other complex systems. Gravitational jostling of orbits by distant moons is comparable to the 
expected radial excursions of small moonlets comparable in size to Bleriot and other propellers in 
the rings. Interactions w ith larger ring p articles can hinder or reverse migration of a larger moon 



on short timescales (e.g.. lTiscarenoll20120 , In all cases, measuring robust migration rates requires 



long time baselines and stable reference frames. 



4.5. Migration in the presence of distant moons: resonances 



In addition to propellers and wakes, Saturn's moons create a variety of structure within the A 
ring. The moons Pan (Encke gap) and Daphnis (Keeler gap) lie within distinct low density gaps 
within the rings. A 7:6 resonance wit h the coorbiting mo ons Janus and Epimetheus is responsible 



for the outer edge of the A ring (e.g.. iLissauer et al.lll985l ). The "ring moons" Atlas, Pandora, and 



Prometheus just outside the ring also produce structure in the ring. 



To illustrate the rich set of other resonances in the A ring (see also ILissauer Cuzzil 119820 , 
Figure [10] shows the result of a simulation of 10,000 tracer particles orbiting Saturn along with the 
major moons (Mimas, Enceladus, Tethys, Dione, Rhea, Titan, and Iapetus) and several smaller 
moons much closer to the ring (Pan, Daphnis, Atlas, Prometheus, Pandora, Janus and Epimetheus). 
Ring particles have large maxima in e around Pan (the Encke gap) and at the outer edge of the 
ring (the 7:6 resonance of Janus and Epimetheus). Other rises in e occur at the orbits of Atlas and 
Prometheus - which lie beyond the outer edge of the ring in regions of very low surface density - 
and at many other labeled resonances. Several of these resonances lie close to but are not obviously 
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coincident with the orbit of Daphnis within the Keeler gap. Curiously, a large ensemble of propeller 
moonlets lie between the 5:4 resonance of Janus/Epimetheus and the 5:3 resonance with Mimas. 

To explore the impact of orbital resonances on migration, we consider whether small moons 
can pass through the 5:4 resonance with Janus and Epimetheus. We set up an ensemble of moons 
with r H = 3-15 km on either side of the resonance near a = 130, 500 km. Because small moons can 
migrate inward or outward, some moons attempt to migrate across the resonance. Others try to 
migrate away from the resonance. In all cases, Janus/Epimetheus stir up the eccentricities of ring 
particles close to the resonance. Large vertical scale height and large eccentricity are barriers to 
fast migration. Thus, we expect the migration of smaller moons to stall close to resonance. 

The results in Figure Qj] confirm our expectations. Resonant interactions with Janus and 
Epimetheus stir up enough ring material to create an obstacle for smaller moons. Moons with 
r H = 3.5 km do not cross the "barrier" set up by the 5:4 resonance. However, moons with r H = 14 
km migrate right through the resonance. 

To understand this resonant barrier to migration, we examine the eccentricity of ring particles 
close to and far away from the resonance. At resonance, Janus and Epimetheus can stir ring 
particles to fairly high eccentricity (e p ~ 10 -5 ; see Fig. [TOj) , Away from resonance, e p is a factor of 
> 10 smaller. At a distance of 130,000 km from Saturn, these ring particles have radial excursions 
of a few kilometers. When they are at resonance, small moons or moonlets with Hill radii smaller 
than a few km are thus embedded in a d ynamically hot ring, which dramatically reduces their 



migration rate (|Bromley Sz Kenyonll2011al . eq. 21). Moons with Hill radii larger than a few km lie 
in a dynamically cold ring whether or not they are within the resonance. These larger satellites 
migrate freely. 

The value of r H that separates these two modes of migration through a resonance depends on 
several factors. With other properties of the rings held fixed, stronger (weaker) resonances produce 
stronger (weaker) eccentricity enhancements, creating barriers for larger (smaller) moons. If the 
strength of the resonance is held fixed, stronger (weaker) damping or viscous spreading weakens 
(strengthens) eccentricity enhancements, enabling (preventing) migration. Thus, understanding 
barriers to migration requires a thorough understanding of the physical processes that impact the 
eccentricity of ring particles. 



4.6. Migration in the presence of distant moons: stirring 



Although distant moons outside of resonance weakly stir the eccentricity and inclination of 
ring particles, t his stirring can have a m ajor impact on the structure of the rings. Using the impulse 
approximation, IWeidenschillingl (|1989l ) showed that distant moons give a particle a small kick in e 
and i every synodic period. For an ensemble of distant moons, he derived a stirring rate 
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where m n , a n and i n are the mass, orbital distance and inclination of the n th moon. These 
expressions are valid for objects on nearly circular orbits. 

For the Saturn system, this approach is valid for moons with separations from the ring, 5a, 
much smaller than the semimajor axis of ring particles. Thus, these equations apply to moons as 
far as Janus and Epimetheus (and possibly Mimas), but are invalid for more distant moons (e.g., 
Titan). With an orbital distance of 1.2 x 10 6 km, Titan tugs more on the barycenter of the entire 
saturnian system instead of stirring the rings. 

Our calculations yield results in rough agreement with the analytic theory. In several simula- 
tions for Figure [TUl the eccentricity of ring particles grows at a rate de 2 /dt ~ 10~ n yr _1 , within a 
factor of two of the prediction from eq. (|37p . With typical e < 10 -6 for ring particles, the growth 
time for stirring by distant moons is a fraction of a year, roughly equivalent to the damping time. 

To compare stirring by distant moons with collisional damping in more detail, we consider the 
energy involved in both mechanisms. Stirring by distant moons results in a kinetic energy input 
per unit mass of V s th ~ de 2 /dt ~2x 10~ 7 dyne g _1 . The power per unit mass dissipated by 
collisions is 

vj 

^coll ~ -rV p TTrlfi oss m p vl (39) 

~ 6 x 1(T 7 dyne/g, (40) 

where m p is the characteristic mass of the ring particles. To derive this estimate, we adopt ring 
particles with r p = 1 m and p p = 1 g cm -3 . The rough equality of these rates suggests that the 
A ring's velocity dispersion and thickness are supported by stirring from distant moons. A precise 
balance between the two energies requires modest changes in the size/density of the ring particles. 

To examine the impact of stirring by distant moons in more detail, we consider simulations of a 
single ring particle interacting with Mimas, the innermost massive moon of Saturn. Each calculation 
begins with a rine particle orbiting at a ~ 136,500 km, with e ~ 3 x 10~ 5 and i ~ 6 x 10~ 5 relative 
to the ring planqj. Mimas begins at a random phase in its orbit, which is tilted by roughly 1° 
relative to the ring plane. 



Figure [T2l shows the results. In a basic set of simulations, Mim as rapidly stirs the e an d i of the 



ring particle (left panels of figure). The growth rate of e follows the I Weidenschillingj (119891 ) estimate 
closely. The growth rate of i is several times faster than predicted. In these simulations, a distant 
perturber (Mimas) continuously forces the alignment of angular momentum vectors between Mimas 
and the ring particle. Thus, the impulse approximation fails and i grows rapidly. Although adding 
collisional damping slows the growth of e and i (right panels), e and i still grow rapidly. 

When we introduce the oblateness of Saturn into these simulations, the growth of e and i 
stalls. In the gravitational field of an oblate Saturn, the ring particle precesses rapidly. Because 
the precession time is short compared to Mimas' orbital period, Mimas cannot force alignment of 
angular momenta. Thus, a single ring particle maintains a constant e and i (Fig. [T2| left panels, 



6 These elements are similar to those of Daphnis. 
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light grey lines). Reducing the precession rate by a factor of 100 enables Mimas some impact on 
the ring particle; e and i then vary periodically (Fig. [T2l left panels, dark grey lines). 

Including damping and oblateness into the calculations produces a different equilibrium be- 
tween stirring and damping. With a reduced precession rate (Fig. 1121 right panels, dark grey lines), 
e and i oscillate about a slowly declining median level. With a normal precession rate around an 
oblate Saturn, e and i decline monotonically (Fig. 1121 right panels, light grey lines). At late times, 
t 3> 10 2 yr, e and i settle into an equilibrium with e ~ i ~ 10~ 7 — 10~ 6 set by the damping rate. 
Larger damping results in a smaller equilibrium e and i. 

This suite of simulations demonstrates that a balance between damping and stirring by distant 
moons sets an equilibrium e and i for ring particles. Throughout the ring, the damping time grows 
in regions of reduced surface density (eq. [8]). Thus, stirring by distant moons is more important 
in regions of low surface density. Because migration slows in regions with large e and i, stirring by 
distant moons can inhibit migration more easily in regions of low surface density. 



5. Daphnis 



If the migration theory outlined in £j2]is correct, the A ring should contain no moon or moonlet 
larger that a few kilometers. Larger satellites migrate to the ring edges on timescales much shorter 
than the ring lifetime. Although some smaller moonlets - like the propellers - can migrate slowly, 
our simulations show that they cannot pass through modest resonances. Aside from Daphnis, the 
ring has no large moonlets or small moons capable of migrating through the resonances shown in 
Figure [TUl Because Daphnis should migrate so swiftly, understanding the origin and the behavior 
of Daphnis is essential to migration theory. 

With an orbi tal distance of 13 6,505 km from Saturn, Daphnis is embedded in the ~ 20 km 
wide Keeler gap (W eiss et al. 2009). It has a modest eccentricity e ~ 3 x 10 -5 and inclination 
i ~ 6.3 x 10 -5 ( Jacobson et al.ll2008l ). Relative to a circular orbit in the ring plane, the radial 
and vertical excu rsions of Daph nis have amplitudes of roughly 10 km, about twice the moon's Hill 
radius of 4.9 km (|Thomaal2010i ). 



Daphnis is a mystery. Despite having a Hill radius within the range identi fied for fast migration , 
its radial drift is ne gligible. Using the measured viscosity of the A ring (e.g., iTiscareno et al.ll2007l ; 
Colwell et al.ll2009i ). the timescale for ring particles to circularize Daphnis' orbit is ^15 00 yr ( eq. [7j; 
see also Hahn 2008). Reducing the inclination to zero follows soon thereafter (see also lHahnll2007l ) . 
allowing fast migration to commence. With an expected total lifetime in the interior of the A ring, 
~ 2500 yr, much shorter than the age of the rings, migration models are closely linked to models 
for the origin of Daphnis. 



5.1. Origins 



There are three possible origins for Daphnis. If this small moon is a recent addition, it might 
be a fragment produced during the impact of a comet or an asteroid with the ring system. Al- 
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ternatively, it might have formed in situ. We reject the idea that the moon results from a three 
(or more) body interaction with Saturn's major moons and outer irregular satellites. The escape 
velocities of Saturn's moons are too low to scatter the moon into the rings. If the current orbit of 
Daphnis is stable for timescales much much longer than the ~ 1000 yr migration time, Daphnis 
is a much older resident of the rings, perhaps a remnant of the tidal event that produced the ring 
system. In this scenario, some set of external processes prevents Daphnis from migrating. To judge 
which of these possibilities is more likely, we consider each in more detail. 



Although there is a low probability that Daphnis is an impact fragment (e.g.. lCuzzi &: Durisen 



1990), this scenario is easy to test with observations. The current orbital eccentricity and inclina- 
tion are appropriate for a small moon settling into the ring (Fig. [6]). The nu merical simulation s 
demonstrate that e and % damp to zero on timescales of 100-1000 yr (see also iHahnl 120071 . 12008). 
For Daphnis' orbit in the A ring, we expect the maximum altitude above the ring to decrease by 
10-100 meters per year. Although there are no current observational limits on this rate, a direct 
measurement would constrain this 'recent external origin' scenario. 



Forming Daphnis in situ is more likely. ICharnoz et al.1 (|2010l . l201ll ) show that many of the 
inner moons of Saturn may have grown from ring material. Although Daphnis lies inside the orbits 
of the moons considered in ICharnoz et al.l (|2010j ) , it is reasonably close to the outer edge of the 
A ring. If Atlas, Prometheus, and other larger moons formed just outside the ring, it is reasonable 
to suppose Daphnis formed with them and migrated inwards to its present location. 

Despite the attraction of this idea, Daphnis' current orbit still requires special timing. The 
size of the Keeler gap surrounding Daphnis' orbit is roughly the correct size for a gap produced 
by a moon like Daphnis. If Daphnis cleared this gap, then it should migrate on short timescales 
(Fig. [5|). Thus, for in situ formation to be viable, it probably happened relatively recently. 



Although this conclusion seems at odds with the ICharnoz et al.1 (120101 ) simulations, observa- 
tions can test this 'recent in situ origin' scenario. Because changes in altitude above the rings and 
the radial motion of Daphnis are detectable, direct measurements can establish whether the orbit 
of Daphnis is damping vertically or migrating radially. 

Because these scenarios seem so unlikely, we consider other processes that might prevent rapid 
migration of Daphnis. As shown in the simulations for Figure [TTT distant moons may establish a res- 
onant barrier which prevents radial migration. In a cursory search for candidates, the Prometheus 
32:31 and Pandora 18:17 resonances seem promising. Although the chaotic behavior of the orbits of 
Saturn's inner moons complicates measuring the strength of resonances and their impact on Daph- 
nis, our simulations demonstrate that resonances can stop the radial migration of moons similar to 
Daphnis. 

Gravitational stirring by distant moons may also prevent migration (Fig. [T2j) . Several of the 
larger moons - notably Mimas and Tethys - have orbits with inclinations i « 1° relative to the 



ring plane. Mimas and Tethys are also locked in a 2:4 mean motion resonance (e.g., I Ailed 1 1969 ; 
Vienne &i Duriez 1992 ; Champenois Viennefl999 ) . Stirring by Mimas and Tethys counter viscous 
damping by the rings. From equations (137p - (138p . we estimate that Mimas doubles the eccentricity 
of Daphnis in less than a decade. With its larger mass, Tethys can double Daphnis' eccentricity 
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in 1.5 yr. Because these timescales are comparable to the damping time, stirring plausibly slows 
down migration. 



5.2. Simulations 

To test the two possibilities to halt Daphnis' migration in more detail, we turn to a set of 
numerical calculations. Currently, we cannot simulate the full set of physical processes which 
impact Daphnis' orbit. Our goal is to show that resonances and stirring by more distant moons 
are plausible mechanisms to prevent migration in Daphnis. 

We begin with a suite of test runs for a single moon with the Hill radius of Daphnis embedded 
withi n the Keeler gap. Starting from Daphnis' current e and i, the moon rapidly aligns with the 



ring (jHahnl 120071 . 120081 ). With a roughly linear decrease in the vertical excursion of Daphnis of 
about 50 m yr -1 , Daphnis' maximum height above the ring plane drops from ~ 10 km to much 
smaller than 1 km in roughly 200 yr (Fig. I13p. Once Daphnis has a vertical amplitude smaller than 
the scale height of the ring, it begins fast migration. 

For the next set of simulations, we consider a Daphnis-sized moon in the 32:31 resonance with 
Prometheus. The resonance pumps up the eccentricity and then maintains a roughly sinusoidal 
variation in e for ~ 10 yr. As e oscillates, the resonance pushes the moon outside the gap. Once the 
moon lies within pristine ring material, damping dominates stirring by the single outer moon. Both 
e and i decline dramatically. When i falls below the threshold if ast < r H /a, the moon commences 
fast migration (Fig. [131 light grey line) . 

Moving Daphnis from the Prometheus 32:31 resonance to the Pandora 18:17 resonance leads 
to similar behavior on longer timescales. Within the Pandora resonance, Daphnis' orbit remains 
inside the Keeler gap for ~ 20 yr. The resonance then moves Daphnis out of the gap into the ring. 
Although it takes awhile for the ring to damp Daphnis' inclination, i eventually drops below ifast- 
Daphnis then begins fast migration. 

When we move Daphnis' orbit to an off-resonance location, its orbit becomes more stable, 
at least on decade- long time scales (Fig. [T3l light grey line). In this set of calculations, Daphnis 
settles into the ring on a time scale comparable to calculations without Prometheus and Pandora, 
or perhaps even faster. There is no clear sign that stirring from these small, nearby moons can 
prevent settling. However the moons do cause some small differences in the evolution of orbital 
elements, which is an encouraging sign: during the first 5000 simulated orbits (7.5 yr), Daphnis' 
inclination shows very little damping. 

These simulations demonstrate the ability of interactions with the distant moons to affect the 
orbital parameters of Daphnis-sized moons on very short timescales. Finding plausible equilibria for 
Daphnis within the Keeler gap requires a detailed parameter search to identify the best combination 
of initial (a, e, i) for the satellite and the outer moons, and (£, f ra d) Td amp ) for the ring. Although 
several test simulations suggest it is possible to find equilibria stable over million year timescales, 
the computational effort is significant and beyond the scope of this paper. 



In a final set of simulations, we consider the reaction of Daphnis to non-resonant interactions 
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with all of the distant moons. The simulations with Mimas in §4.5 show that a balance between 
damping and stirring from a single moon can maintain a large inclination, with % > if as t- Stirring 
by Prometheus and Pandora can reduce settling rates, at least over several thousand orbits. Thus, 
stirring from all the outer moons could maintain a large inclination, probably with some jitter due 
to the constantly changing gravitational field along Daphnis' orbit. 

As in Figure [TOl these simulations include all of the major moons and the ring moons. The 
coordinates of the m oons and Daphnis relative to Saturn come from the JPL Horizons systems 



(jGiorgini et all 11999 ) 00:00:00 GMT, 2012-June-22; we adjusted the ring plane to enable starting 
conditions for Daphnis similar to those in Figure [T3l The calculations include the A-ring with an 
empty Keeler gap, damping and diffusion of ring particles, and Saturn's oblateness. 

Figure [14] shows the results of one simulation. As in Figure [T3l the thin black curve shows 
the decay of e and i of a lone Daphnis orbiting within the Keeler gap and interacting only with 
ring particles. Orbital interactions with the full ensemble of Saturn's moons yields relatively rapid, 
small amplitude fluctuations in a and e. This jitter is similar in amplitude to the jitter of any small 
moon in the ring (see Fig. [9]). On much longer timescales of a few months to half a year, Saturn's 
moons produce 0.5 km oscillations in the z max of the moon. On decade timescales, <z max >, the 
median height of the moon above the ring plane is roughly constant. Thus, stirring by Saturn's 
moons appears capable of maintaining the large inclination of Daphnis relative to the A ring. 

The predicted oscillations of Daphnis are sensitive to several characteristics of our calculations. 
In our approach to ring dynamics, ring particles have no self gravity and therefore do not react to 
changes in local surface density. Furthermore, the ring plane itself is fixed and does not evolve in 
response to torques. Despite these limitations, the calculations make a robust prediction: as the 
mass of a small moon grows and as the surface density of ring particles around the satellite declines, 
the distant moons of Saturn produce an oscillatory behavior in the height of the small moon above 
the ring plane. These oscillations can prevent small moons like Daphnis from migrating through 
the ring. 

The predicted quasi-periodic oscillations in the altitude of Daphnis above the ring plane are 
detectable. The predicted amplitude of 200-400 m is a significant fraction of the typical maximum 
altitude of roughly 9 km. The quasi-period of a few months yields a rate z max ~ 0.4-0.8 km yr -1 



a few times larger than the rates of motion derived for the propeller moon Bleriot (jTiscareno et al 



201u ) . Thus, accurate observations of Daphnis' inclination or maximum height above the ring plane 



would test our interpretation. 



5.3. Summary 

Our results suggest several plausible origins for Daphnis. A 'young' Daphnis can result from an 
impact fragment or from growth beyond the outer edge of the A ring. In both of these recent origin 
scenarios, ring material on either side of the Keeler gap should damp Daphnis' orbit on short, ~ 
1000 yr, timescales. Interactions between Daphnis and the rings should also induce fast migration. 

An 'old' Daphnis has been resident in the rings for timescales exceeding ~ 10,000 yr. Because 
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external processes prevent Daphnis from migrating, Daphnis' residence time in the ring could 
approach the ring lifetime of > 0.1-1 Gyr. Our simulations show that orbital resonances with 
distant moons can impact Daphnis' orbit. Identifying the combination of physical properties for 
Daphnis, the rings, and the distant moons needed to maintain Daphnis in its current orbit requires 
a comprehensive suite of simulations that is beyond the scope of this paper. 

Our simulations also demonstrate that orbital interactions between Daphnis and distant moons 
(e.g., Mimas and Tethys) outside of any resonance can impact Daphnis' orbit. Maintaining Daphnis 
in its current orbit requires a fine balance between stirring by the distant moon(s) and damping by 
ring particles outside the Keeler gap. 



5.4. Observational Tests 

Observations can test all of these scenarios. If Daphnis is young, we predict detectable, mono- 
tonic changes in the orbital parameters (a, e, i). If Daphnis is old, changes in (a, e, i) depend on the 
importance of resonances. Outside resonance, we expect modest changes in phase with the motions 
of distant moons such as Mimas. Within resonance, the behavior depends on the strength of the 
resonance relative to damping by the ring. Although measuring orbital parameters in a chaotic 
reference frame is challenging, direct measurement of Daphnis' altitude above the ring plane and 
position in the Keeler gap provide equivalent tests of our predictions. 



6. Conclusions 



Saturn's A rin g is rich in str ucture, with an amazing, complex variety of transient and long- 
lived features (e.g.. lEspositoll2010l ). The smallest ring particle s lie in a matrix of se lf-gravity wakes, 
which form coherent features on scales smaller than 100 m (IHedman et al.l 120071 ). The 150 or so 
moonlets embedded in t he ring have physical radii between 50 m and 250 m, and are located in 
three ~ 1000 km annuli (jTiscareno et al.ll2008l ). The small moons Pan and Daphnis orbit within 
distinct gaps with very low surface density. Together, these moons, moonlets , and smaller particles 
provide an excellent testbed for migration theory (see also ICrida et al.l l2010i : iPan &: Chiang! 120101 : 
Rein & Paoaloizoul I2OI0I : ISalmon fc Charnozll201o1 : lTiscarenol[20rj . and references therein). 



As outlined in Figures Q] and [H theory predicts several distinct modes of migration through 
the rings. In the very viscous A ring, moonlets with r H < 0.1 km migrate on timescales much 
longer than the age of the solar system. Large moons with r H > 20-30 km also migrate very slowly. 
Between these limits, migration timescales are shorter than the ring lifetime. For r H ~ 2-20 km, 
the timescales are comparable to a human lifetime. 

In our simulations of Saturn's A ring, fast migration of small moons with r H ~ 2-20 km is 
ubiquitous. When a small moon lies within a gap, it eventually pulls enough material from the gap 
edges to begin fast migration (Figs. [3H4J) . When a small moon is embedded within a continuous 
ring, it eventually clears enough material from its corotation zone to begin fast migration (Fig. [5]). 



In a pristine ring of small particles, numerical simulations show that moonlets and small moons 
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migrate at the rates predicted by analytic theory (Figs. [5]andE]). In numerical simulations that 
include damping, gravitational interactions between small moons and ring particles, and viscous 
spreading, the moons clear their corotation zones on timescales T c i ear . When the viscous timescale 
is longer than T c i ear , these moons open a gap in the ring and begin to migrate. As long as r H < r^t, 
migration is fast. Otherwise, migration is slow. 

In Saturn's rings, damping and viscous spreading smooth out the wakes left behind by mi- 
grating moons. Thus, multiple moons can easily migrate rapidly through the rings. This situation 
contrasts with protoplane tary disks, where stirring b y migrating protoplanets overcomes damping 
by smaller planetesimals ([Bromley Kenvonll2011bl ). 



Several physical processes within the rings limit or halt migration. Interactions between 
small moons and (i) very large ring particles or (ii) local enhancements in ring surface den- 
sity lead to stochastic migration, often reducing migration rates significantly (jCrida et al.l 120101 : 
Rein &: Papaloizoull2010l : lTiscarenoll2012l ). Migrating moonlets and very small moons cannot pass 
the 'resonant barriers' formed by moons outside the A ring (§4.4). Pairs of resonances can 'corral' 
these objects and prevent them from migrating in or out through the ring (Fig. [TO]) . Stirring by 
distant moons outside resonance can also overcome damping and prevent migration (§4.5). 



These results lead to several clear conclusions regarding migration in the A ring. 



Prop eller moonlet migration. With Hill radii less than r. 



ap 



2 km (Tiscareno et al 



20081 ). the propeller moonlets should experience slow, embedded migration with radial drift 



2012), and oscillatory motion from interactions with corotating material (jPan &: Chiang 



201 







rates below 10 m yr 1 ( Fig. [2|). Stochastic migration in a self-gravitating ring ( Rein fc Papal oizou 
201C ; ICrida et al.ll2010i ). interactions with surface density variations in the ring ( Tiscareno 



2012 ) have similar low drift rates. For most moonlets, measurements of radial drift limit the 
rate to below 0.1 km y r -1 . At least one larger moonlet, Bleriot, shows non-Keplerian motion 
(ITiscareno et al.l 120101 ) , with a drift rate that varies in time. Our new results suggest that 
the radial drift of Bleriot and other moonlets can be confined by orbital resonances. Indeed, 
most propellers appear to be confined to the 'corral' between the Mimas 5:3 resonance and 
the Janus/Epimetheus 5:4 resonance (Fig. 1101) . 



The propelle r moonlet size dist r ibuti on. Moonlets have a range of sizes up to about 
r H ~ 1.5 km (jTiscareno et al.l 120081 . |2010j). Migration theory predicts an upper limit of 



r gap ~ 2 km. Larger objects migrate out of the ring on timescales of 100-1000 yr. The 
good agreement between r gap and the maximum size of propellers is a strong confirmation of 
migration theory. 

Non-Keplerian motion of moonlets. Moonlets orbit Saturn in a gravitational field which 
is constantly modified by much more massive moons orbiting at larger distance. The constant 
jostling of moonlets by distant moons results in chaotic radial variations in distance from 
Saturn (Fig. [9]). Aside from complicating the measurement of a moonlet's orbital elements 
within our simulations, this motion may complicate interpretat ions of non-Keplerian motion 
of moonlets within the ring (e.g., Bleriot, ITiscareno et al.ll2010i ). 
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The 'young' Daphnis. If some external physical mechanism cannot maintain Daphnis in 
its current orbit, it must be a recent, ~ 10 3 — 10 4 yr, addition to the ring. Migration theory 
predicts rapid damping in (e, i) followed by fast migration. Instruments on board Cassini can 
measure the predicted damping in height above the ring plane. 

The 'old' Daphnis. If orbital resonances or stirring by distant moons can maintain Daphnis 
in its current orbit, it may be very old. Prom our simulations, we identify sources of stirring 
which might maintain Daphnis' orbit for timescales exceeding 10 4 yr. Measuring oscillations 
in (a, e, i) can place strong constraints on this mechanism. 

Pan. With r H = 19 km at an orbital distance of 133,584 km, this sm all moon orbits within the 



300 km-wide Encke gap (jCuzzi k, Scargldll985l ; |Jacobson et al.ll2008l ). Pan's Hill radius is very 



close to rf as t. As with Daphnis, orbital resonances and stirring from distant moons probably 
prevent any form of migration. In addition, these effects may operate to help maintain the 
Encke gap, whose edges are relatively far, ~ 8r H , from Pan. The torques from ring material 
in this configuration are too weak to cause measurable radial drift. 

• Atlas. Orbiting within the Roche Division at the outer edge of the A ring, Atlas' Hill 
radius, r H = 22 km, is also close to rf ast . Although it is too far away from ring material to 
have any chance at migrating, Atlas may have migrated out through the ring to its present 
location. As with Pan and Daphnis, understanding the internal physics of ring material and 
the interactions with distant moons might yield clues about its origin and recent history. 

• Pan and Atlas. For both moons, detecting or limiting variations in (a, e, i) can provide 
important constraints on the physical mechanisms within the ring and interactions between 
ring particles, these small moons, and distant moons. 

To conclude, our analytic estimates and numerical simulations demonstrate that Saturn's rings 
are an important laboratory for testing migration theory. Current results are encouraging. Within 
the A ring, we can explain aspects of the radial and size distributions of propellers and we iden- 
tify plausible mechanisms for limiting migration in Daphnis. We also identify several important 
observational tests which can be accomplished with current satellites. 

Aside from including additional constraints on ring properties as new observations provide 
them, future studies of ring dynamics will eventually require includin g the ring self-gravity directly. 



Self-gravity adds to the rich dynamical structure of the rings (e.g., iLewis Stewart! l2009l ). which 



likely has interesting consequences for the migration of moons and moonlets (e.g., Fig. [9|). 
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Appendix 



~ 4r H of a small 



The onset of fast migration relies on the presence of some ring material within 
moon. For a moon within a gap, fast migration proceeds only when the half-width of the gap is 
< 4r H . Because the relationship between the moon and the size of the gap it creates is critical to 
migration, we consider several approaches to assess this relationship. 

U), the viscous torque on ring material 



In the torque-balance approach (e.g.. iLissauer et al. 



1981 



as it diffuses into a gap of widt h Aa is set equal to the to rque at the gap edge applied by the satellite. 
From weak-scattering theory (|Lin &: Papaloizoul Il979l ) , this balance yields a simple relationship 
between the gap half-width and the satellite's Hill radius: 



Aa 



GM\ 1/6 
a 5 u 2 J 



(1) 



If the time to fill the gap is the viscous time (e.g., Tgn in eq. [7] in the main text), it is possible to 
eliminate the viscosity in this expression and derive a minimum gap size, as in eq. ([2]) for r gap . 

An alternative approach equates the radial distance ring material can diffuse in a synodic 
period, yVT syn , with the displacement in semimajor axis a particle receives after every close en- 
counter with the satellite (proportional to the right-hand-side of eq. ([1]) above; see also eq. [5] of 
Bromley & Kenyon 2011b). This balanced is achieved at an orbital distance from the satellite 



Aa 



(GM_\ 
V J 



1/18 



A/3 



(2) 



This result is similar to the result in eq. (pQ). 



More detailed treatments derive some form of the radial diffusion equation acr oss the gap and 
solve explicitly for t he surface density distribution and the half-width of the gap (e.g. jHourigan &; Ward 



1984 



Rafikov 



2001 



Although these solutions yield much better estimates for the surface density 
distribution, results for the half-width of the gap differ from the simpler treatments by less than a 
factor of two. Given the uncertainties, either eq. ([1]) or eq. ([2]) provides a reasonable approximation 
to the half-width. 

To apply these two results to Saturn's A ring, we adopt the nominal parameters in eqs. (H 
(1161) from the main text and derive 



Aa ~ 0.97 



for balancing the torques and 



Aa ~ 1.4 
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(4) 



for balancing scattering and viscous diffusion. For small moons capable of fast migration, r H ~ 
2-5 km, these approaches yield similar results for Aa. For much larger moons, torque balance 
consistently yields a larger half-width. 
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For the rest of this section, we focus on half-widths derived by balancing scattering and viscous 
diffusion. Requiring Aa > 4r H yields the maximum Hill radius for a satellite capable of drawing 
material from the gap edges, ~ 23 km, when scattering balances viscous diffusion. This limit is 
comparable to rf ast . 

In principle, moons with r H > 23 km can prevent material from encroaching within 4r H . These 
moons migrate slowly in the type II mode. Because gravitational scattering falls off so steeply with 
orbital distance from the satellite, viscous diffusion constantly tries to drive material into the gap. 
If diffusion produces a higher surface density and viscosity at the edges of a gap, the gap will shrink, 
enabling the satellite to draw material away from the gap, possibly initiating fast migration. 

To compare the simple theory with Daphnis in the Keeler Gap and Pan in the Encke Gap, we 
adopt r H ~ 5 km for Daphnis and r H ~ 20 km for Pan. For Daphnis, the balance between scattering 
and viscous diffusion occu rs at Aa ~ 12 km , somewhat smaller than the observed half width, ~ 13- 
20 km, of the Keeler gap (jWeiss et al.ll2009l ). For Pan, the balance is at Aa ~ 75 km, much smaller 
than the observed half-width, ~ 160 km, of the Encke Gap. Based on this analysis, we identify 
Daphnis as a good candidate for fast migration. Pan is a poor candidate for fast migration, at least 
under current conditions in the Encke Gap. However, barring other effects such as resonances with 
distant moons, Pan may be unable to halt long-term diffusion of ring material to within 4r H . 

In our calculations, large moons achieve a balance between viscous diffusion and gravitational 
scattering. The numerical results support the limits derived in eq. (|4|). Moons with r H > 20-25 km 
can migrate in the type II mode. Smaller moons with r H ~ 2-20 km may undergo fast migration. 



REFERENCES 

Alexander, R. D., & Armitage, P. J. 2009, ApJ, 704, 989 
Allen, R. R. 1969, AJ, 74, 497 

Artymowicz, P. 2004, Debris Disks and the Formation of Planets, 324, 39 
Borderies, N., Goldreich, P., & Tremaine, S. 1984, ApJ, 284, 429 
Bromley, B., & Kenyon, S. J. 2006, AJ, 131, 2737 
Bromley, B., & Kenyon, S. J., 2011, ApJ, 731, 101 
Bromley, B., & Kenyon, S. J., 2011, ApJ, 735, 29 

Bryden, G., Chen, X., Lin, D. N. C, Nelson, R. P., & Papaloizou, J. C. B. 1999, ApJ, 514, 344 

Champenois, S., & Vienne, A. 1999, Celestial Mechanics and Dynamical Astronomy, 74, 111 

Champenois, S., & Vienne, A. 1999, Icarus, 140, 106 

Cook, A. F., & Franklin, F. A. 1964, AJ, 69, 173 

Charnoz, S., Salmon, J., & Crida, A. 2010, Nature, 465, 752 



- 30 - 

Charnoz, S., Crida, A., Castillo- Rogez, J. C, et al. 2011, Icarus, 216, 535 

Colwell, J. E., Nicholson, P. D., Tiscareno, M. S., et al. 2009, Saturn from Cassini-Huygens, 375 
Crida, A., Morbidelli, A., & Masset, F. 2006, Icarus, 181, 587 

Crida, A., Papaloizou, J. C. B., Rein, H., Charnoz, S., & Salmon, J. 2010, AJ, 140, 944 

Cuzzi, J. N., Lissauer, J. J., & Shu, F. H. 1981, Nature, 292, 703 

Cuzzi, J., Clark, R., Filacchione, C, et al. 2009, Saturn from Cassini-Huygens, 459 

Cuzzi, J. N., Burns, J. A., Charnoz, S., et al. 2010, Science, 327, 1470 

Cuzzi, J. N., & Durisen, R. H. 1990, Icarus, 84, 467 

Cuzzi, J. N., & Scargle, J. D. 1985, ApJ, 292, 276 

D'Angelo, C, Dapnis, R. H., & Lissauer, J. J. 2010. larXiv:1006.5486l 

D'Angelo, C, Kley, W., & Henning, T. 2003, ApJ, 586, 540 

D'Angelo, C, & Marzari, F. 2012, ApJ, 757, 50 

Esposito, L. W. 2010, Annual Review of Earth and Planetary Sciences, 38, 383 
Gladman, B. 1993, Icarus, 106, 247 

Giorgini, J. D., Yeomans, D. K., Chamberlin, A. B., et al. 1996, Bulletin of the American Astro- 
nomical Society, 28, 1158 

Goldreich, P., & Tremaine, S. D. 1978, Icarus, 34, 227 

Goldreich, P., & Tremaine, S. D. 1978, Icarus, 34, 240 

Goldreich, P., k Tremaine, S. 1979, ApJ, 233, 857 

Goldreich, P., & Tremaine, S. 1980, ApJ, 241, 425 

Goldreich, P., & Tremaine, S. 1981, ApJ, 243, 1062 

Goldreich, P., & Tremaine, S. 1982, ARA&A, 20, 249 

Hahn, J. M. 2007, ApJ, 665, 856 

Hahn, J. M. 2008, ApJ, 680, 1569 

Hahn, J. M., & Malhotra, R. 1999, AJ, 117, 3041 

Hahn, J. M., Spitale, J. N., & Porco, C. C. 2009, ApJ, 699, 686 

Hedman, M. M., Nicholson, P. D., Salo, H., et al. 2007, AJ, 133, 2624 

Hourigan, K., & Ward, W. R. 1984, Icarus, 60, 29 



Ida, S., Bryden, G., Lin, D. N. C, & Tanaka, H. 2000, ApJ, 534, 428 
Jacobson, R. A., Antreasian, P. G., Bordi, J. J., et al. 2006, AJ, 132, 2520 
Jacobson, R. A., Spitale, J., Porco, C. C, et al. 2008, AJ, 135, 261 
Kenyon, S. J., & Bromley, B. C. 2004, AJ, 127, 513 
Kenyon, S. J., & Bromley, B. C. 2008, ApJS, 179, 451 
Kenyon, S. J., & Bromley, B. C. 2009, ApJ, 690, L140 

Kirsh, D. R., Duncan, M., Brasser, R., & Levison, H. F. 2009, Icarus, 199, 197 

Levison, H. F., Morbidelli, A., Gomes, R., & Backman, D. 2007, Protostars and Planets V, 669 

Lewis, M. C., & Stewart, G. R. 2000, AJ, 120, 3295 

Lewis, M. C., & Stewart, G. R. 2009, Icarus, 199, 387 

Lewis, M., Stewart, G., Leezer, J., & West, A. 2011, Icarus, 213, 201 

Lin, D. N. C., & Papaloizou, J. 1979, MNRAS, 186, 799 

Lin, D. N. C., & Papaloizou, J. 1986, ApJ, 309, 846 

Lissauer, J. J., & Cuzzi, J. N. 1982, AJ, 87, 1051 

Lissauer, J. J., Shu, F. H., & Cuzzi, J. N. 1981, Nature, 292, 707 

Lissauer, J. J., Fabrycky, D. C., Ford, E. B., et al. 2011, Nature, 470, 53 

Lissauer, J. J., Goldreich, P., & Tremaine, S. 1985, Icarus, 64, 425 

Lubow, S. H., & Ida, S. 2010. larXiv: 1004.41371 

Malhotra, R. 1993, Nature, 365, 819 

Masset, F. S., & Papaloizou, J. C. B. 2003, ApJ, 588, 494 
Masset, F., & Snellgrove, M. 2001, MNRAS, 320, L55 
Morbidelli, A., & Crida, A. 2007, Icarus, 191, 158 

Morbidelli, A., Tsiganis, K., Crida, A., Levison, H. F., & Gomes, R. 2007, AJ, 134, 1790 
Ohtsuki, K. 1992, Icarus, 98, 20 

Paardekooper, S.-J., Baruteau, C, Crida, A., & Kley, W. 2010, MNRAS, 401, 1950 
Paardekooper, S.-J., Baruteau, C, & Kley, W. 2011, MNRAS, 410, 293 
Pan, M., & Chiang, E. 2010, ApJ, 722, L178 
Pan, M., & Chiang, E. 2012, AJ, 143, 9 



-32- 

Pan, M., Rein, H., Chiang, E., & Evans, S. N. 2012, MNRAS, in press (|arXiv:1206.3583j) 

Papaloizou, J. C. B., Nelson, R. P., Kley, W., Masset, F. S., & Artymowicz, P. 2007, Protostars 
and Planets V, 655 

Perrine, R. P., Richardson, D. C, & Scheeres, D. J. 2011, Icarus, 212, 719 

Porco, C. C, and the Cassini Imaging Team, 2005, IAU Circ, 8524, 1 

Porco, C. C, Weiss, J. W., Richardson, D. C, et al. 2008, AJ, 136, 2172 

Pringle, J. E. 1981, ARA&A, 19, 137 

Rafikov, R. R. 2001, AJ, 122, 2713 

Rein, H., & Papaloizou, J. C. B. 2010, A&A, 524, A22 

Salmon, J., & Charnoz, S. 2010, Bulletin of the American Astronomical Society, 42, 981 

Salmon, J., Charnoz, S., Crida, A., & Brahic, A. 2010, Icarus, 209, 771 

Salo, H., & Karjalainen, R. 2003, Icarus, 164, 428 

Salo, H., & Schmidt, J. 2010, Icarus, 206, 390 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Showalter, M. R. 1991, Nature, 351, 709 

Showalter, M. R., Cuzzi, J. N., Marouf, E. A., & Esposito, L. W. 1986, Icarus, 66, 297 

Shu, F. H., Cuzzi, J. N., & Lissauer, J. J. 1983, Icarus, 53, 185 

Syer, D., & Clarke, C. J. 1995, MNRAS, 277, 758 

Tanaka, H., Takeuchi, T., & Ward, W. R. 2002, ApJ, 565, 1257 

Thomas, P. C. 2010, Icarus, 208, 395 

Tiscareno, M. S., Burns, J. A., Nicholson, P. D., Hedman, M. M., & Porco, C. C. 2007, Icarus, 189, 
14 

Tiscareno, M. S., Burns, J. A., Hedman, M. M., et al. 2006, Nature, 440, 648 

Tiscareno, M. S., Burns, J. A., Hedman, M. M., & Porco, C. C. 2008, AJ, 135, 1083 

Tiscareno, M. S., Burns, J. A., Sremcevic, M., et al. 2010, ApJ, 718, L92 

Tiscareno, M. S. 2012, Planetary and Space Science, in press (arXiv:1206.4942). 

Tiscareno, M. S. 2013, in Planets, Stars, and Stellar Systems, Volume 3: Solar and Stellar Systems 
Oswalt TD, French L, and Kalas P, eds. (Dordrecht: Springer), in press (arXiv:1112.3305). 

Vienne, A., & Duriez, L. 1992, A&A, 257, 331 



- 33 - 

Ward, W. R. 1988, Icarus, 73, 330 

Ward, W. R. 1997, Icarus, 126, 261 

Ward, W. R., & Hourigan, K. 1989, ApJ, 347, 490 

Weidenschilling, S. J. 1989, Icarus, 80, 179 

Weiss, J. W., Porco, C. C, k Tiscareno, M. S. 2009, AJ, 138, 272 
Youdin, A. N., & Kenyon, S. J. 2012. larXiv:1206.0738l 



This preprint was prepared with the AAS IATgX macros v5.2. 



-34- 




Fig. 1. — A qualitative comparison of predicted migration rates as a function of Hill radius for 
satellites in a particle disk. The solid curves show rates in viscous disks; the dashed curves apply- 
in the limit of zero viscosity where r gap — the Hill radius of a satellite that can clear a gap — is 
set by the disk scale height, h. A satellite drifts radially according to the larger of the fast and 
embedded migration rates for its Hill radius. The dip in the embedded migration rate for f ra( j > 
illustrates that migration may be suppressed in a viscous disk if r H is in the "transgap" regime, 
where the satellite is larger than the disk scale height but not large enough to open a gap. The 
curve representing fast migration shows an upper plateau with a migration rate independent of 
r H , a steep decline above the limiting Hill radius rf ast (see text), and a similarly sharp attenuation 
in da/dt below r H = h, where the typical disk particle eccentricity e rms becomes important. The 
shaded region is suggestive of radial drifts from stochastic migration within some fixed time period. 
If the source of the stochasticity is clumping of disk particles on scales of A c i ump , then the effective 
rate falls off as 1/V H . Otherwise the effective rate is constant. 
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Fig. 2. — A schematic of predicted migration rates as a function of Hill radius for satellites in 
Saturn's A ring. The solid curve shows the upper limit of the predicted rate for a wide range of 
satellite masses. The dashed line is the low-viscosity limit of embedded migration. For r H > ffastj 
the solid line below the dashed line reflects the viscosity-dominated Type II mode. The shaded 
region represents stochastic migration, bounded by an approximate average yearly radial drift 
ass uming that the rin g con tains clumps of 30 m (thes e rates are drawn from a broad range reported 
by ICrida et al.1 (|2010l ) and iRein &; Papaloizoul (|2010l )). The data points and shaded region for the 
propeller moonlets indicate Hill radii only; with the exception of the moonlet Bleriot, migration 
has not been detected. 
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Fig. 3. — Fast migration of small moons in a ring of massive particles with S = 40 g/cm 2 , a total 
annular width of 56 km, and a gap with a halfwidth of 14 km, which covers the corotation zone. 
Simulations begin with the moon precisely centered in the gap. Each curve represents the change 
in orbital distance, Aa, of a moon with r H = 7 km and an initial semimajor axis of 130, 000 km. In 
each simulation, the particles in the ring have identical mass, set to 10 -2 (lightest shading), 10 -3 , 
10 -4 , 10 -5 , or 10~ 6 (darkest shading) times the mass of the moon. Solid (dashed) lines indicate 
inward (outward) migration. Although migration is more stochastic with larger ring particles, the 
migration rate is nearly independent of the masses of ring particles. When a moon reaches the edge 
of the ring, migration ceases. 
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Fig. 4. — Snapshots of small moon migration in a closing gap. The top image showa the location 
of a Daphnis-size moon (r H = 4.9 km, the central white disk of 2.3 km in radius, the approximate 
physical size of Daphnis) at Daphnis's orbital distance (136,505 km) when it is initially on a circular 
orbit in a pristine ring with uniform surface density S = 40 g/cm 2 , centered on a clear annulus the 
size of the Keeler gap (full width of 40 km). The ring particles (masses equal to 0.1% times that of 
the satellite) are indicated with a brightness that correlates with density. The lower three images 
show views from the same orbital distance from Saturn but at times equal to 46.6 yr (second from 
top), 48.2 yr, and 49 yr. Radial diffusion (z^ ra d = 85 cm 2 /s) causes 'the ring to slowly creep into 
the gap, but once particles get close, fast migration initiates and the moon heads toward Saturn. 
Each image represents an area of 800 km in width and 120 km in height, oriented so that Saturn is 
located toward the top of the page. The "exposure time" for each image is several orbital periods, 
to enhance the structures in the ring. 
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Fig. 5. — The disk mass in the corotation zone and the radial drift of a small moon with r H = 14 km. 
The starting orbital distance is 130,000 km. The upper plot shows the instantaneous disk mass m co 
inside an annulus of half- width 2r H ; the lower plot is the orbital distance of the moon relative to 
its starting position. The black curve refers to calculations for the moon and ring alone; the light 
curve indicates results for calculations which include Saturn's prominent moons (Fig. I10|) . In both 
instances, migration stops after ~ 170 years, when the moon hits the inner edge of the simulated 
disk. 
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Fig. 6. — Orbital settling of captured small moons. Simulations of moons with Hill radii of 10 km 
(light gray curves), 14 km (dark gray curves), and 18 km (black curves) show the evolution of 
orbital elements in an initially uniform disk with A ring parameters. In all cases, the time for each 
moon to become embedded in the disk is very short. 
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Fig. 7. — The onset of fast migration of a captured moon. The traces show the change in semimajor 
axis (upper plot), eccentricity (middle plot) and inclination (lower plot) of a small moon with 
r H = 6.5 km in a disk with £ = 20 g/cm 3 . The initial orbit of the moon is the same as for 
Daphnis; the disk properties are modified for computational feasibility. The disk has an 80 km- 
wide gap (similar to the Keeler gap along Daphnis' orbit in terms of the moon's Hill radii). The 
disk particles' eccentricity and inclination damp as expected for the A ring. The viscosity, however, 
is turned up by a factor of 15, thus reducing the viscous spreading time in the gap by a factor 
of 3 compared to the Keeler gap. The simulation shows that the eccentricity damps quickly, but 
migration does not begin until the inclination is also small. After migration starts, its rate is 
comparable to the theoretical prediction (dashed curve) , and it halts once the moon nears the edge 
of the simulated disk. 
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Fig. 8. — Fast migration of small moons in a ring of massive particles with £ = 30 g/cm 2 . Each 
curve represents the change in orbital distance, Aa, of a moon with r H = 14 km that initially has a 
semimajor axis of 130, 000 km. In each simulation the particles in the ring have identical mass, and 
are set to either 10, 50, 100, 300 and 600 times smaller than the mass of the moon. The darkest 
curve shows the case with the most extreme mass raio (600:1); lighter shades have progressively 
less extreme mass ratios. The lightest curves, corresponding to runs where the particle masses are 
10% of the moon mass, show stochastic behavior. 
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Fig. 9. — Fast migration of a small moonlet in the presence of Saturn's massive moons. In these 
simulations a moonlet with r H = 0.7 km starts off at an orbital distance of 130,000 km in a ring 
with surface density S = 0.1 g/cm 3 . The black curve corresponds to ring particles that are 1% 
of the moonlet mass, while the dark gray curve shows results for ring particles that are twice as 
massive. The corotation zone is initially clear, and no viscous diffusion can fill in the gap. Thus, 
fast migration ensues, and there is a subsequent change in orbital distance, Aa (solid black curve; 
the dashed line shows the theoretical prediction). Note that the true orbital distance (not shown) 
is given by the sum of Aa and the orbital fluctuations that result from orbital "jostling" of the 
moonlet by the distant moons. We estimate this jostling using a moonlet with identical starting 
conditions as the other simulations, but without any ring particles (light gray curve). 
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Fig. 10. — The eccentricity of simulated ring particles as a function of semimajor axis in the 
presence of massive moons. The simulation started from cold initial conditions and evolved tracer 
particles over a 4-year period. Some of the strongest resonances are labeled; off-resonance pumping 
of eccentricity is also apparent throughout the ring. The black circles designate moons; the size 
of each circle i s suggestive of the rel ative mass. The open circles are a select sample of propeller 
moonlets from lTiscareno et al.l (|2008). None of these symbols are intended to quantify eccentricity, 
just the semimajor axis. The gray bands are the zones where propeller moonlets seem to reside 
exclusively. The figure also labels the location of the Encke and Keeler gaps, as well as the inner 
edge of the Roche Division (containing the light gray points). 
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Fig. 11. — Fast migration in the presence of resonances from distant moons Janus and Epimetheus. 
The two panels show the output of simulations with isolated small moons on circular orbits at 
a position near the resonances. The black-colored traces show the radial position of each moon 
Aa = km relative to an orbital distance of a = 130, 500 km. For reference, the light grey lines give 
the location of the resonances from Janus and Epimetheus (upper and lower lines, respectively). 
The five runs in the left panel show fast migration of a smaller (r H = 3.5 km) moon; the runs in 
the right panel are for a larger moon (r H = 14 km). For each set of runs the initial conditions 
are identical, except that the moon and ring are radially offset by a differing amount. The results 
indicate that the smaller moons cannot cross the resonance; the larger ones can pass through. 
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Fig. 12. — Stirring of a single massless ring particle by a distant moon (Mimas). The ring particle 
has initial orbital elements similar to Daphnis. Left panels: stirring by Mimas without collisional 
damping. Right panels: stirring by Mimas with collisional damping. In each panel, the black curves 
show the growth in e and % when Saturn and Mimas have a spherically symmetric gravitational 
potential. Stirring grows continuously. When Saturn's oblateness is included in the calculations, 
the ring particle precesses rapidly. Stirring and precession yield a rough equilibrium in e and i as 
expected for a ring particle. When precession is slowed by a factor of 100, the variations in e and 
i become periodic due to the precession of the ring particle's angular momentum vector precesses 
(dark gray curves). 
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Fig. 13. — Simulations of Daphnis' orbit: damping, stirring and resonances. The traces show 
changes in the semimajor axis (upper panel), eccentricity (middle panel) and maximum altitude 
above the ring plane (lower panel, the inclination is i ~ z max /a) of a simulated moon with r H = 
2.6 km at a = 136, 50 5 km. The initial eccentricity and inclination are similar to those of Daphnis 
(jJacobson et al.ll2008l ). To help isolate physical effects, the oblateness of Saturn is not included in 
these simulations (see Fig. [T2"j) . The curves illustrate the evolution of (i) a small moon embedded 
in the Keeler gap (black curve), (ii) a small moon with its position adjusted to coincide with the 
32:31 resonance of Prometheus (light gray curve), and (hi) a small moon, off-resonance but stirred 
by Prometheus and Pandora (dark gray curve). Without resonance or interactions with distant 
moons, Daphnis settles into the ring (~ 10 2 — 10 3 yr) and then begins fast migration. In a single 
resonance, Daphnis is scattered out of the Keeler gap and begins to migrate. Evolution within the 
Pandora 18:17 resonance is similar but on a longer timescale. Off resonances, the orbit is stable on 
timescales of at least 10 4 yr before Daphnis begins to settle into the ring. 



-47- 



(0 



N*. 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



Daphnis only 

Major moons, ring moons 








4 6 
time (yr) 



8 



Fig. 14. — Simulated evolution of Daphnis' orbit when stirred by Saturns moons. The moon orbits 
within a gap in surface density with a density contrast and width similar to the Keeler gap. As 
in Figure [T3l the black curve shows changes in the moon's orbital elements solely from interaction 
with the A-ring; over the course of 6000 orbital periods, Daphnis slowly settles in the ring and 
eventually begins to migrate rapidly. The light curve shows the evolution when Daphnis interacts 
with ring particles a nd all of the moons o f Saturn. Using initial phase-space positions from JPL's 
Horizons calculator (|Giorgini et al.lll996l ). Saturn's moons produce quasi-periodic fluctuations in 
the inclination, which prevent Daphnis from settling into the ring plane. Thus, interactions with 
distant moons plausibly prevent fast migration of Daphnis through the A ring. 



